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NASA Technical Memorandum 78238 


Computer Program Development and User's Manual for Program PARACH 


I. INTRODUCTION 

This user's manual furnishes complete documentation on Program 
PARACH, a digital simulation currently being used to investigate the 
interaction dynamics of a parachute-spent Solid Rocket Booater (SRB,) 
during terminal descent. Program PARACH is written in Fortran and is 
operational on the Univac 1108 and .EMR 6050 computers. The FMR 6050 is 
used because of the graphic display capabilities of the computer. 

Section II presents a discussion of the mathematical model and a 
summary of resulting equations, while Section III furnishes a description 
of the program's subroutines. Section IV is a listing of input variables 
used in the program and their definitions. Finally, Section V gives a 
program flow chart and a listing of the program. 


II. MATHEMATICAL MODEL 

The math model simulates a SRB xjith two parachutes attached at 
different points on the SRB, Runs made up to now only use one parachute. 
The other one is zeroed out. 

Assumptions - Several assumptions are made to simplify the 3 
dimensional math model formulation and its implementation into a digital 
computer simulation called Program PAEJ,JH. These are; 

„ All bodies are considered as rigid bodies, 

. The mass and aerodynamic forces of the suspension lines 
are neglected, 

. The center of pressure and center of volume for the 

parachute canopy are assumed to be located at the center 
of mass of the canopy material. 


Parachute Configuration - The program regards the parachute con- 
figuration as a rigid body. It has six degrees of freedom; three 
translational and three rotational. The main disturbance is due to 
aerodynamic forces and moments. Use of chute apparent mass and inertia 
properly distinguishes the parachute from a conventional rigid body. 

Booster Configuration - This program also regards the booster as 
a rigid body. It has six degrees of freedom, 3 translational and 3 
rotational. There are aerodynamic forces and moments computed on this 
rigid body. 

System Configuration - The system configuration considers the three 
rigid bodies to be joined by frictionless attachments that transmit no 
moments. The attach point constraint relationships are formed by requiring 
the inertial velocities of the bodies to be the same at the attach point. 

Atmospheric Mass Density - Density is calculated using the following 
equation: p = K * EXP (Aq + Aj^ Z + A 2 + A^ + A^ Z^ + A 5 Z^). 

The capability to use a constant density and emit the equation is also 
available. 

Parachute and Payload Drag - Drag is a function of velocity, mass 
density of air, angle-of-attack , reference area, and other factors. Wake 
effects on either body were not considered. 

Rotational Damping Coefficient - The damping on payload and 
parachute due to rotational velocities is currently set to zero. How- 
ever, the program includes the provisions to include the rotational 
damping derivatives. 

Retro Thrusters on the Payload - The program has the capability of 
having eight retro-rockets on the payload, positioned and skewed as 
desired. There can be more if the appropriate dimension statements are 
enlarged. At ignition, each thruster follows the thrust profile in the 
THR table. 

Wind Model 

Table Look-Up (Subroutine TBL) - A wind table is used as an 
alternative to a programmed wind profile. The wind table uses sub- 
routine TBL which is an interpolation scheme. The subroutine is used 
in every other table look-up in this program. 
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Parachute Side Force - Side forces are considered in this 
program as part of the parachute dynamics. It is presented here as a 
function of the angle-of-attack (see figure 1). 


Program Equations - The following program equations are 
computed, accounting for inertial, reaction, gravity, aerodynamic, and 
thruster forces. The program equations are computed for both the booster 
and the parachute yielding both translational and rotational accelerations 
and the sum of the forces acting on the two bodies. 


Q, = + F^^^ + F^^^ + F^^^ + 


Where: - an 18 x 18 diagonal matrix whose elements 

are the mass and inertia values for the SRB and 
the two parachutes, 

F^^) = Inertial Forces and Moments (18 x 1 vector) 

F' = Reaction Forces and Moments (18 x 1 vector) 
F^^^ = Aerodynamic Forces and Moments (18 x 1 vector) 
= Control Forces and Moments (18 x 1 vector) 

= Gravity Forces and Moments (18 x 1 vector) 


cf = (Ul, Vi, Mj, Pi, Qi, Rj. Uj. ''2. Wj, P2, Qa. ^ 2 ’ “3. ''3. “3. '’ 3 > ^ 3 : 


^ 3 > 


The following equations are used to compute the Inertial 
Forces and Moments (F^^)), 
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and: Uj^ , Vj^ , are the linear velocities of the 

respective bodies, is the mass term and are the angular 

velocities of the respective bodies. 

Coordinate Transformations: 

The general transformation between any two 
Euclidean coordinate systems is written in terms of Euler parameters 
(quaternions ) as : 
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G. 

a. 


(q[)^ - (q^)^ - <qp^ + (qj)' 
2 qj - q3 qj) 

2 (qj q3^ ^ q^ qj) 


tl- /-|i _L 


n / i i i i\ 

2 (q2 qs - q^^ q4> 


2 (<ll <13 - “2 

2 (q2 q3 + qj’ qj) 

- (qj;)^ - (q^)^ + (q^)^ + 



Where: 

G]^ (transformation from inertial to SRB frame) 

G2 (transformation from inertial to parachute 
A frame) 

G3 (Transformation from inertial to parachute 
B frame) 

G? = G2 g'^ (transformation from SRB to parachute 

o 

G^ = G3 G^ (transformation from SRB to parachute 

B) 

The q's in the coordinate transformation matrices 
are quaternions given by the following equations : 

i CO il) 0, . li) .. 9. (D. 

Oi = sin ^i cos ^i cos 1 - sin ^1 sin 1 cos ^1 


A Q ih 

= sin cos ^ 

2 2 


iii . ip. , . ib . . <p. e. 

^1 cos ^1 + sin “1 sin 1 cos 1 


. 


q^ = sin ^ cos ^ cos ^ - sin _i sin ^ cos ^ 

^ 2 2 2 2. 2 2 

i lb Q (P . lb .. Q .. (P . 

q, = cos ^ cos _i cos ^ + sin sin ^ sin ^ 

2 2 2 2 2 2 


Where 0^, cp^, and i/)^ are the Euler angles 
i = 1 refers to the SRB 
i = 2 refers to parachute A 
i = 3 refers to parachute B 


5 


•P' H' 



The differential equation governing the evo- 
lution of the quarternions as a function of time is given by: 
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Likewise, the Gravity Forces and Moments 
are derived from the following: 
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Where: F^^^ = 
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m^ is the appropriate actual mass term 
g is the gravitational force 


The thrust vector due to retro-rockets is 
computed using the following equations 
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f 

= THR(t) . 2 DCOSX(i) 

= THR(t) ^ DGOSY(i) 

i^= 1 

= THR(t) 2. DCOSZ(i) 

= THR(t) X ( YR(i) DCOSZ(i) - 

= THR(t) J ( ZR(i) DCOSX(i) 

= THR(t) ^ ( XR(i) DCOSY(i) - 

i = 1 

= 0 


ZR(i) DCQSY(i)) 
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Where: THR(t) is the thrust profile as a 

function of time (t). 

DCOSX(i), DCOSY(i), and DCOSZ(i) 
are the direction cosines of the 
thrust vector. 

XR(i), YR(i), ZR(i) are the moment 
. arms of the retro-rockets. 

The aerodynamic forces and moment vectors are 
computed using the following equations: 
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The calculation of the reaction force is quite 
a bit more complicated as far as computational procedures are concerned 
and is given by the following equations 
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FI6.1. PROGRAM COORDINATE SYSTEMS 



III. PROGRAM SUBROUTINES 


Subroutine ; RUNGE 

Purpose ; Subroutine RUNGE is a subroutine used to integrate the 

differential equations and produces a table of integrated 
values , 

Use ; PiUIGE (KUTTA, TIME, DT, NVAR, NDVAR) 

Description of Parameters ; 

KUTTA - is a control integer (controls the number of times through 
integration loop) 

TIME - is time of integration 

PT - is the time increment 

NVAR - is the number of variables to be integrated once 
NDVAR - is the number of variables to be integrated twice 
Remarks : Fourth-order RUNGE-KUTTA 

Subroutine ; TBL 

Purpose ; This subroutine provides a linear interpolation between data 
points in a tabular function of one variable. 

Use : TBL (X, Y, XI, N, Yl) 

Description of Parameters ; 

X - is the independent variable 

Y - is the dependent variable 

XI - is the X argument 

N - is the number of data points in the table 

Y - is the result desired 


Remarks : The X array of data must be in ascending order. 
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Subroutine; LINSYS 


Purpose ; This subroutine takes the inverse of a matrix. 

Use ; LINSYS (A, NROWSA, NCOLSA, DET, B, ARRAY) 

Description of Parameters ; 

A - is the array containing the coefficients of the linear system. 

NROWSA - is the number of rovs of the square Matrix A. 

NCOLSB - is the number of columns of the Matrix B in the equation 
A*X = B. 

DET - is the determinate of Matrix A 
B - is the Matrix B = A*X 

NARRAY - is the number of rows of the arrays in which A and B are 
stored in the calling program. 

Subroutine ; XT AN 2 

Purpose ; This subroutine is used to define correct quadrant. 

Use ; XTAN 2 (A, B) 

Description of Parameters ; 

A - is the value of the numerator of the angle. 

B - is the value of the denominator of the angle. 

X 
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IV. INPUT PARAMETERS 


F0RTRAN 

MNEMONIC 

ENGINEERING 

SYMBOLS 

TYPE 

DEFINITIONS 

XI 


R 

Inertial Location of SRB CG in 
the X Direction 

Y1 


R 

Inertial Location of SRB CG in 
the y Direction 

Z1 


R 

Inertial Location of SRB CG in 
the z Direction 

U1 


R 

Inertial Velocity of the SRB CG 
in the x Direction 

VI 

Vl 

R 

Inertial Velocity of the SRB CG 
in the y Direction 

W1 


R 

Inertial Velocity of the SRB CG 
in the z Direction 

X2 

X2 

R 

Inertial Location of Parachute A 
CG in the x Direction 

Y2 

^2 

R 

Inertial Location of Parachute A 
CG in the y Direction 

Z2 

^2 

R 

Inertial Location of Parachute A 
CG in the z Direction 

U2 

CM 

R 

Inertial Velocity of Parachute A 
CG in the x Direction 

V2 

V2 

R 

Inertial Velocity of Parachute A 
CG in the y Direction 

W2 

W 2 

R 

Inertial Velocity of Parachute A 
CG in the z Direction 

Note; 

Any consistent set 

of units 

may be used. 
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FORTRAN 

MNEMONIC 

ENGINEERING 

SYMBOLS 

TYPE 

DEFINITIONS 

X3 

X 3 

R 

Inertial Location of Parachute B 
CG in the x Direction 

Y3 

^3 

R 

Inertial Location of Parachute B 
CG in the y Direction 

Z3 

Z3 

R 

Inertial Location of Parachute B 
CG in the z Direction 

N 

N 

I 

(9) Number of Data Points in Table 

N1 

Nl 

I 

(41) Number of Data Feints in Table 

U3 

U3 

R 

Inertial Velocity of Parachute B 
CG in the x Direction 

V3 

^3 

R 

Inertial Velocity of Parachute B 
CG in the y Direction 

W3 


R 

Inertial Velocity of Parachute B 
CG in the z Direction 

PI 


R 

Angular velocity of SRB About the 
xl Axis 

Q 1 

Ql 

R 

Angular velocity of SRB About the 
yl Axis 

Rl 

Rl 

R 

Angular velocity of SRB About the 
zl Axis 

THIl 

ex 

R 

2nd Euler Angle of SRB 

PHIl 

.h 

R 

3rd Euler Angle of SRB 

PS 11 

ii)i 

R 

1st Euler Angle of SRB 

P2 

P2 

R 

Angular Rotation of Parachute A 
about the x2 Axis 
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F0RTRAN 

ENGINEERING 


MNEMONIC 

SYMBOLS 

TYPE 

Q2 

^2 

R 

R2 

R 2 

R 

THI2 

©2 

R 

PHI2 

02 

R 

PSI2 

02 

R 

P3 

^3 

R 

Q3 

Q3 

R 

R3 

^3 

R 

THI3 

O 3 

R 

PHI3 

03 

R 

PSI3 

03 

R 

XIXl 

^xi 

R 

XIYl 

^yi 

R 

XIZl 


R 

XIX2 

^^2 

R 


DEFINITIONS 

Angular volocity Parachute A 
About the y2 Axis 

Angular velocity of Parachute A 
About the z2 Axis 

2nd Euler Angle of Parachute A 

3rd Euler Angle of Parachute A 

1st Euler Angle of Parachute A 

Angular velocity of Parachute B 
About the x3 Axis 

Angular velocity of Parachute B 
About the y3 Axis 

Angular velocity of Parachute B 
About the z3 Axis 

2nd Euler Angle of Parachute B 

3rd Euler Angle of Parachute B 

1st Euler Angle of Parachute B 

Moment of Inertia About XfAxis 
of SRB 

Moment of Inertia About y^Axis 
of SRB 

Moment of Inertia About Zj^Axis 
of SRB 

Moment of Inertia About X 2 Axis 
of Parachute A 
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F0RTRAN 

ENGINEERING 


MNEMONIC 

SYMBOLS 

TYPE 

XIY2 


R 

XIZ2 

^2 

R 

XIX3 


R 

XIY3 

^y3 

R 

XIZ3 

^23 

R 

SI 

Si 

R 

S2 

S 2 

R 

S3 

S 3 

R 

YL 

1 

R 

YPL 

2 

R 

YPL3 

3 

R 

XXW2 

Mig 

R 

XXW3 

M2g 

R 

XMl 

Ml 

R 

XM21 


R 

XM22 

^2 

R 

XM23 

Mz 2 

R 

* Mass is 

equal to total 

mass plus 


DEFINITIONS 

Moment of Inertia About y 2 Axis 
of Parachute A 

Moment of Inertia About Z 2 Axis 
of Parachute A 

Moment of Inertia About X 3 Axis 
of Parachute B 

Moment of Inertia About Axis 
of Parachute B 

Moment of Inertia About ^ Axis 
of Parachute B 

Surface Area of SRB 

Surface Area of Parachute A 

Surface Area of Parachute B 

Reference Diameter of SRB 

Reference Diameter of Parachute A 

Reference Diameter of Parachute B 

Weight of Parachute A 

Weight of Parachute B 

*Mass of SRB 

*Mass of Parachute A in x Direction 
*Mass of Parachute A in y Direction 
*Mass of Parachute A in z Direction 
apparent mass 
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FORTRAN 

ENGINEERING 



MNEMONIC 

SYMBOLS 

TYPE 

DEFINITIONS 

XM31 


R 

Mass of Parachute B in x Direction 

XM32 

•V3 

R 

Mass of Parachute B in y Direction 

XM33 


R 

Mass of Parachute B in z Direction 

G 

g 

R 

Gravity 

XLXl 


R 

xl Distance From CG of SRB to 
Attach Point of Parachute A 

XLYl 

Si 

R 

yl Distance From CG of SRB to 
Attach Point of Parachute A 

XLZl 

Lzl 

R 

zl Distance From CG of SRB to 
Attach Point of Parachute A 

XLX2 

S2 

R 

xl Distance From CG of SRB to 
Attach Point of Parachute B 

XLY2 

CM 

R 

yl Distance From CG of SRB to 
Attach Point of Parachute B 

XLZ2 

S2 

R 

zl Distance From CG of SRB to 
Attach Point of Parachute B 

XLAXl 

Sa 

R 

x2 Distance From Attach Point 
to CG of Parachute A 

XLAYl 

^yA 

R 

y2 Distance From Attach Point 
to CG of Parachute A 

XLAZl 

^zA 

R 

z2 Distance From Attach Point 
to CG of Parachute A 

XLBXl 

LxB 

R 

x3 Distance From Attach Point 
to CG of Parachute B 
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F0RTRAN ENGINEERING 


MNEMONIC 

SYMBOLS 

TYPE 

DEFINITIONS 

XLBXl 

^xB 

R 

y3 Distance From Attach Point 
to CG of Parachute B 

XLBYl 

^yB 

R 

z3 Distance From Attach Point 
to CG of Parachute B 

XLBZl 

^zB 

R 

Direction Cosines of the Eight 
Retro-Rockets (I = 1,8) 

DCOSX(I) 

D cos x(i) 

R 

Direction Cosines of the Eight 
Retro-Rockets (I = 1,8) 

DCOSY(I) 

D cos y(i) 

R 

Direction Cosines of the Eight 
Retro-Rockets (I = 1,8) 

DCOSZ(I) 

D cos z(i) 

R 

Distance From CG of SRB to Retro- 
Rockets 

LCX(I) 

x(i) 

R 

Distance From CG of SRB to Retro- 
Rockets 

LCY(I) 

y(i) 

R 

Distance From CG of SRB to Retro- 
Rockets 

TIMX(I) 

t 

R 

Time Table 

THR(I) 

"^hr 

R 

Thrust Table 

Z(I) 

Z 

R 

‘Altitude Table 

UWD(I) 

Uv 

R 

u Wind Table 

WD(I) 

Vw 

R 

V Wind Table 

AL(I) 

CL 

R 

Alpha Table 

FCNAl (I) 


R 

Table of Aerodynamic Forces 
Acting or SRB (Normal) 
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F0RTRAN 

MNEMONIC 

ENGINEERING 

SYMBOLS 

TYPE 

DEFINITIONS 

FCAAL(I) 


R 

Table of Aerodynamic Forces 
Acting or SRB (Axial) 

FCMAl(I) 


R 

Table of Aerodynamic Forces 
Acting or SRB (Moment) 

FCNA2(I) 

W 2 

R 

Table of Aerodynamic Forces 
Acting on Parachute A (Normal) 

FCAA2(I) 

^Aa2 

R 

Table of Aerodynamic Forces 
Acting on Parachute A (Axial) 

FCMA2(I) 

C%2 

R 

Table of Aerodynamic Forces 
Acting on Parachute A (Moment) 

FCNA3(I) 

C%3 

R 

Table of Aerodynamic Forces 
Acting on Parachute B (Normal) 

FCAA3(I) 


R 

Table of Aerodynamic Forces 
Acting on Parachute B (Axial) 

FCMA3(I) 

^%3 

R 

Table of Aerodynamic Forces 
Acting on Parachute B (Moment) 
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V. FLOW CHART AND PROGRAM LISTING 
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TRANSFORMATION 
OF VECTOR (fi) INTO 
INERTIAL FRAME 
AND PRINT 



INTEGRATE THE 
INERTIAL VELOCITIES 
AND PRINT 



Graphics Display for "PARACH" 
(Sample Results) 


The following plots were made on the ERR-6050 computer. 

These are included to illustrate the graphics display capability 
of the computer. 

The curves shown are representative results for one parachute 
with the appropriate parachute characteristics listed at the top 
of each graph. 
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S C bCOlsS- THE NUHBER OF COLUHNS OF THE HATRIX B IW THE EQUATION 

f C A*X : B. 

7 C . — That is, TMlS.ROUlJMLJil Bl USED. TO SOI WE THE. ABOWC EWUAT IO»L 

a C WHEN B <aNO XI are either SInGlE VEC|0RS or SEwERAL vectors 

9 C ..conjoined IN.A HATRlIi,. NOIE THg^T |e. J .iS A U»,l T JlATfilRj eHC 

10 C SOLUTION X IS THE INVERSE OF THE MATRIX A. 


H. 

1? 

13 

lA 

IS 

It 

.11 

If 

19 

1C 

21 

22 

<3 

2* 

15 

26 

27 

IF 

29 

3T 

31 

32 

33 
3A 
3S 
26 
37 
3F 
39 
AC 
Al. 
A2 
A3 
AA 
A5 
A6 
A? 
Af 

A9 

SC 

51 

52 

53 

5A 

55 

56 


pet is the DE TIRHIWANT OF THE HATRIX A . PET SHOULD BE CHECREO TO 


c 

c 

SEE IF IT Is ZERO, FOR IN THIS CASE THE HATRIX IS SINGULAR AND THE 
RESULTS RETURNED ilBC MEANINGLESS, ExCERT FOR.OET* 


c 

c 

6 IS THC HATRIX WHICH CONTAINS ThE RIGHT HAND SlOL OF THE ABOVE 
.EQUATION. IT IS ALSO. THE LOCaTION OF THE SOLUTION hATRIXi.A, ON 


t 

c 

EXIT FROM THE SUBROUTINE. 

NAhRAY IS the number OF ROWS OF THE ARRAYS IN WHICm A AND B ARE 


c 

c . 

STORED IN THE CALLING PROGRAM. THIS HAY BE EQUAL TO NROiiSa, 
BUT NEED NOT BE. THUS THE USER HAT STORE A ( j | JN. A L ARGER .. 


c 

c 

arra' Than is strictly necessary* this is useful when the exact 

DIPENSION OF. A. IS ilNUNOWN 4.T £OnR 1L*T10N tlME*. 



INTEGER Z |, Z2,Z 3, 7^,25, 26, 2 T, 2B, 29, Z 10, Z1 1,212,213,219 
DIMENSION Alil.PIlT 



N = NROWSA 
NN r NCOLSB 



NK = NARRAY 
DET=1 , 

IJIOOOSO 


DO 2 K:l,N 
2l:K«|K-l |*NK 

1J10U060 

Tu 

IFIK-N) 30,3,30 
CONTINUE 



TCST:ABS lAtZill 
KPirR.l 


L=K 

— fiOAiraPiiJi 


22=1*IK-| I9NR 

1FITEST>ABS 1 A | 22 II 1 3 1 , 9 . A 


n 

TES1=ABS IAI22II 


A 

CONTINUE 
IFTL-AI 91,3,91 


Al 

DO S J : N,N 
23=L*IJ-1 l*NK 



TEHP:AI2',I 
29ZR*|J-1 I9N9 


. 5_. 

AI23I:a(29I 
AI29 IzTEhP 



DO IS J=1,NN 
25:L*IJ-1 l•NI« 



TEMPrplZSI 

26ZK.IJ-1I9NK 


IS 

bl2SI=B(26| 

BI26I=tENP 


3 

OET r -OET 
DLT=0ET*A 1211 



A(2II=1./A(21I 


REPRODUCBILri'Y OF TH 
ORIGINAL PAGE IS POOR 






r- 


ilMSTS 0*TC 0«I03?^ 


57 

58 

59 
6D 
bl 

„fc2_ 

IMK-NI 35,19, 3S 

35__ _ DO 6 J = NPl.N _ 


Z7rN« 1 J>1 I9NK 
A(Z7I-A iZ7i9«iZ11 


6 Continue 

19 nO 16 1 = 1 . NM 


b3 

Z9=a* IJ-l I9N9 


69 

1*. BIZ9I=BIZ9>9«IZ1I 


65 

00 7 1:1, N 


66 

Zlu:i*U-ll9Na 


67 

im-Mi 160,7,160 


bf 

160 rAC=«uioi 


69 

00 a j=kpi,n 


7tl 

Z11=1*IJ-1I9NN 


71 

Z12=K*IJ-I »*N8 


17 

_ irAlZllI 


73 

0=C>FAC*A|Z12I 


,19. 

1M1.0E-9*ABS lEl-ABS 1 0 1 1 162. 16 f . 1 6 l 



TS 161 0 = C< 


16 U,l »U11I=0... 

IT 8 continue 

le . oo.ie. j .=. I ,NN 

79 Z13=I«( J-t l*NK 

hC ___iil=IL±ljTll£JUS 

at 18 eizi3izBizi3i-r«c*Biziai 

6? 7 continue 

S3 2 CONTINUE 

aa _ return 
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iIN7Lf,*TPFS(bl.HAlN 

_ BE*l. _NPAR» 

REAL KAPP ,KOV| ,K0V2«K0V3iH2D0T 

.RE<kL^iCHl 

real LCX|8|,LCViei,LCZ(8l 

real mi 8j 18 i,MH3.3l.H2t3.3I.H3l3.3).HHl8.1ftl.mNVt6.6».KU>,»> HOOOQOOC 

dimension PFR7I5U) ,PFR8I50| ,PrR9|50» ,PFRL7I 1 2 I , PFRL 8 1 1 ^ » ,PF RL9 U 2 » 

_ l,PS2(SLI.PS2LU2l«^l>Oti&All&U)tPQBAlL(12I.PQeAR2(&Oj»PQeA2LII2.Li _ 

2P22CSO) ,PZ2LU2i tPOXHSOi.PDXlLl 12) tPALl ISO>tPALIL(I2) • 

3PRi«5o» tPRU«l2) tPPT1150)iEPTIUl2liP*L2«5Q)tPAL2J.J_I2», _ 

APOl (501 ,P0IL 1121 ,PD2HSol fPDZlLi 121 1 

iP PH50I. PPlUi2) 

dimension PtIME ISO I tP2l ISO) •TIMCL 1 12),2iL( 12) 

.. tl^’ENSlON 2CEA21i6)>2CNA2<36)t2CAA2«36)tPERAL2U2)4«AZiJ2J.jNjli2J2J 

9t2XA7l2) 

. DIMENSION ZCNAUsI* 2CHMA5 » tZC aaI ^ 5 ) » XMA AL 1 18 1 • i»A2.12 > ,2X A I Z1 ,NA2 « 2 


DIMENSION BBI3, 3),X0MI 18, 1 ) , XOM U ( 3 1 1 ) , X0N2 1 1 3 , 1 ) 

. DIMENSION 6113.34^iU4-K}i»DlU»}JtD1211i?U(lljUL2liE3l3il) 

dimension XL1I3,3),xL2(3,3) ,XLA1I3,3| XLB1|3,3) ,B11(3,3),B|2(3,3) 

DIMENSION Bl3>3.3) ,BlA13 t3),BlS< 3j, 3) ,6 16 1 3, 3 ) ,B21 1 3 j 3j {622 i3 , 3 ) 

DIMENSION B23I3,3),B2AI3,3),B2SI3,3),B26(3,3),XOM1I3,3),xOM2(3,3) 

D IJi£N 5 1 ON XOM 313.3).BD1H3.3).0D1213.3 ).rD1313.31.BD1M3,3) 

DIMENSION B0I5|3,3),8D16I3,3),B021(3,3),B022I3«3)>B023I3,3) 
dimension B021l3t3l,B0Z5JJj3),B028l5,3),BIC,18),B016,^l8LtAA)3t3) 
DIMENSION X0M3K 3, 1) ,X«n 3,3),Xy2<3*3),xv3(3,3),Xll|3,3),X12(3,3) 
DIMENSION XI3U.3) , F1H3.1 ) ,F 2H 3, 1 ) . F 3H 3, 1 ) , XL 1 H 3 , 1 ) , xt2 H 3 , 1 ) 
DIMENSION XL3I(3,l),FII18,l),XGll3,i),XG2l3,i),X63(S,i),F6i«3,l) 
DIMENSION FG2t3.1).FG3l3.1),FGn8.l).V l M13il)tV2M«3il>,»JM13,ll 
DIMENSION FAI18,l),7IAl),WD|Al),UlNll3,l),UIN2<3,l),UIN3<3tll 

dimension At 115) iFCNA2ll S>tFCNA3ll5).FCAA2(15).FCAA3llS). 

lFCMfc2ll5),FCMA3(15) 

DIMENSION X0MD118 .1).FS Hietl).FR<18,ll 

DIMENSION AiAi6,i8),BIBI 18,6),FSi 18,1) ,FR1I 18,r),FR2il8,r) 

DIMENSION X OMT > 18. 18l .XOMFt 18. 1) 

DIMENSION 0C0SX(8),DC0SY(8),0C0S2I8),FC( 18, 1 ) , T INX I 9 ) , THR I 9 ) 

COMMON/FUNT/V<39)jOyMli^) 

COMMON/FIRT/Fi 39 > •0UM2«21 ) 

EQUIVALENCE 12XA2»l )iPER) tl2XA212),AL2) 

EQUIVALENCE 1 2X A 11 ) , M ACH 1 ) , 1 2X A i 2 > • AlI ) 

EQU IVALENC E IXQMDl H , F t 1 ) ) ■ 1 00 1 1 . F « 19) ) ,100 1 2,F) 2 0 ) ) , I QOl 3 , F < 2 1 ) ) 
ll0021,Ff22)),IQD22,F(23)),(0023,FI2A)),lQD31,Fi2S)),l0D32tFI26)), 
_?I0033,FI27) ),IQDAljFI28))ilODA2jFI29)) ,I00A3,FI30)) 

EQUIVALENCE ( 0 1 1 , V 11 9 ) ) , 1 0 12 , V I 20 > ) , ( Ql 3 , V < 2 1 ) ) , 

1 1021 ,T<22 ) ) , «Q2?ill21iJj*023tV I2A) ) , ( 0 31 , V( 25 ) ) , I 032jt V I 26 1 ) , 
2I033.YI27 I ) , |QA1,YI28 ) ) ,Ioh 2,V«29) ) ,IOn3f TI3U) ) 

equivalence <Ym.xOMll).U l), |XOM<2) .Vl),IXOM 13 1,Ml),IXOM«H),XOMl 

*ll,l),Pl),IX0ril5),Cl),(X0M(6),Rl),IX0M(7),'u2|,(X0M(8),V2),(X0M(9), 
_ •M2) , I XOMI ID) tX0M21 1 Ij 1 LtP2»jlX0'’J II • »02) , » XOM U2X,R^2 jiiXOM 11 3 ) ,1^3 
• , IXOMI lA ) ,V3) , ixOMlis) f M3), IXOMI 1 6 ) , XOM 3 1 U , 1 ) , P3 ) , I X0 m117 ) ,qS) , 

7lxlO,F I 31 1 ), (YlOfF J32) )j I210,F 1 33) ) , 

7IX20,Fi 3 A ) ) , Iy2D,F I 35) ) , (22P,F I 36) ) , 

7 1X3D,F 137 ) ) ,1 Y 3D.F 1 38) ) , 1Z3D ,F 139) ) , 


MOOOygSo 

M0000060 

M0000070 

M0000080. 

M000009a 

MQQDOlOO 

MOOOOllO 

M0000120. 

MOOOOl3o 

MOOOOIAO 

MOOOOISO 

M0000160 

M0000170 


M0000200 

hOUOOZIo 

M0D0 O22D 

M0000230 

M00002A0' 


,M00002S0 

MQ000260 

M00002 7U 

M0d00280 

M000Q290 

Hboou3od 

1 M0000310 

,MUqoo320 

)M00pU330 

MOOOQ3AO 


7IX1,Y131)),)Y1,V13Z)),1Z1,V133)), 

7IX2,YljA)),lV2,Y13S))tlZ2,V136)), 

7|X3,Vl37)),lV3,vi38)),lZ3,Yi39)),IXOMtie),R3) 



57 equivalence (FA 1 17, 1 1 ,VN3> , IFA I 16 . 

58 llFAl l.D.FxD.lFAi 2 ,1 1 ,FY 1) , IFAI 

1> 

3 

,YN3), 

ltlL.F2, 

»,lFAl_Atl»jTU». 

S9 IIFAI 5, 1 1 ,YMl > , IFA 1 6 ,' 1 1 • yN 1 1 , 1 F A 1 

1 

’.ll.FxJ 

),|FA1 6,1),FY2), 


6 b 

61 

.62 

63 


_ UIF*I 9.U .r22l i!EUlUxJJ-iVL2lt irAUl,}i,yH2).|fJli2tUxIN2)4 
llFA(l3tl>>FX3lt(F*<16tl>tFV3ltlF«US.l>»FZ3l.tF*U6,l).YL3l. 

iiFG i j ,j J ,EG1 Ii .in.irei z.ii.fGu z.m.i FGt 3.HifsH3.n». 

llFGt T,U,F62|1 In.lFGI 8 , 11 , FG2( 2 * 1 1 I . I FG ( 9 , 1 1 , F 62 I 3 1 1) I f 

J.(FG(13.1>iFG3llil)it.i£611*ijl^tf.63|2^1) l«lFGtlL*ii-iLe.}iJfil*-i 

IIFII l,ll,Flllltin,(Fll 2.1I.FlII2,in.(Fll 3, 1 | ,F1 1 1 3 ,1 1 I , 

.. HFIl 7>li.E21(ifILUi£li-9tniF2l(2ill)f<Fil ?.lKF^It3.1|l, 

l(Flll3,ll,F31(i,in,(FIUii,ll,F3II2tlll,(Flll5.1i|F3l(3,n>t 

HFH 6.ii,xLiii i.iii.iFH s.ii .xliHz.d » i«f i« 6 u>t?.kU13i^j-Lt- 

llFniOfllfxL2I(l.in.lFl,ll,l|,XL2I(2,ll),IFI(12,i),XL2ll3tlll« 

_ HFHJ6.i»,XL3iUiiUiJiiiniUixL3i«2imnFi(i8,i)jXL3n3,i»i 

LOU I valence (viNii,ii,un.iviHi2tiifVii.iviHi3.n.Hii 

LOUlVALENCt (V 2FU1 • U tU2i i ( V2M (2. 1 I « V2 > t ( V2N) 3 • M tW2L 

equivalence IV3H(l»lltU3),IV3Hf?,ll,V3lt(V3HI3tl)tW3l 
.LQulV 

,_2 tX0K31(2tlttfi3)>IX0N31l3tniR3) 

DATA PFRL7/6HPaI> Rltll*6H / 

_ DATA PFRL8/faHPAP R2t.ll66H / _ 

DATA PFRL9/6HPAR R3,H*6H / 

DAT A P S2L/6 HS 2 . n«6H i 

DATA POBAlL/GHOeARl ,11*6H / 

DATA pQBA2L/6HQ8AR2_tli*6H / 

UAT« PZ2L/6HZ2 .11*6H / 

DATA PPlL/bHPl_ tU* 6H / 

DATA P01L/6HQ1 tll*6H / 

DATA PR1L/6HR1 .U*6H 


DATA PDX1L/6H0X1 ,ll«6N / 

DATA POY1L/6HOV1 tli* 6H _/ 

Data pozil/6h6zi •ii*6h 7 


DATA 

PAL1L/6HAL1 

,1166H 

/ 

data' 

PAL2L/6HAL? 

,11*6H 

/ 


data ;1L/6H71 tH*6H 

DATA !_THR_m ,I=l, 9 l^^ji 3 _ 9 V 0 D.i 39 V 00 . , 0 . , 0 .i 0 ., 0 .,g.jU,/ 

DATA ( Al I i I .i= 1 , 15 ) /O.. ib. tZO. t 30 . (VO. tSO. f 60 . t 70 . t 80 > t 90 . t 

100 . ( 1 1 0 , 1 1 20 * 1 1 30 3 ^ 0 . / 

OATAlOCOSXI I) ,I : l, 8 l/ 86 u./ 


DATA tOCOSZI I ) ,T = 1.8I/8*-!./ 

UA7A «LC* 

DATA ILCV (I),l=i.8l/8*.b/ 

pA!*_«LC2.llJ,lllj8l/8*.C/ _ _ 

DATA IFCAA2II),l:l,15»7.55,.52S,.*i7S,.3eS,.27S,0.,0.,0.fU.,0..0. 

»-t P j, ±0»z 

DATA |FCNA2III,I=1 ,15)/0...0G, .125t.l9,.25,10*.3/ 

DATA «FCMA2ll)il;l(15)/g.,-.g25,-.g625t-.095,-.125,10*-.15/ 

data (ZCAA(n,i:i,15>7.7S,I.2S,l775,2.35,2.35,2.0,i.5,U.8S,0.22~, 

*-.25,r.25,.5i.22,-_.15tr._L5/_ _ _ 

DATA (ZCNAIII,iri,iS)/0.0,.5,1.5,3.i,S.4,7.6,B.2,6.A,6.,6.,6.,6. 


LATA (ZCHA|l)flrXfl5)/O.Uffl«8f3.,3.3,6.7,6.7t8.C,6.0,1.7f*2.S, 

6*5.8 ,*8.2 ,*10. ,*10 .,."5. / ____ 

DATA IZCAAm,iri6,3oi7l. 35, 1 . 5 , 1 . 9 , 2 . 5 , 2 . V S , Zl 0 ,1 . 5 , . 9 , . 3 5 . • . 1 , 


0«TE D«ti379 


page 


2 


MQ 003 e^ 

N 00 b 0390 

HPUDOVOO 

NooboAib 



»4«iH 



DATE 090379 RAGE 3 

IIS 

♦•*3f •25f*#l f*»75/ 

U*TA UCN^lll Lf 6,3.7 ^6.4.8 



116 

117 

*Ii.0iiu.Mf9.e.a.**f 9*0/ 

DAT* IZChAIII *1 = 16. SOI/O.Oi 1*6.2. S. 3. 0.3*6, 6. 8. 8.2.9. 0 1.2.8 i>7.Qi 



llfl 

IIS 

•-I0.7.-12.9.-13.0.-11.8.-9.7/ 

DAT* 12 CNAI 11*1 = 31 *951/0*. *7. 2*. 9*1. 7*5, 10. 7. 13*. 19*25, 15*8. 16*. 



120 

121 

915*6, 19. 7, 19.4.11*5,9 *5/ 

DATA |2CMAllJ,l= 31 .951 /Q , . 1 .8* ?. 3*3* 9. 5« Li 7 *_} .7*8.5. 3.U* .-6. .*11 * 7 



H2 

123 

12*1 

J25 

i', -17 *',-20 *7, -19 *1,-8*/ 

DAT* t7C3AHl,l=31f951/l*25.1*7,2*,2*2.2**l*65.1*05..25.-.b5i-*6t- 



♦•5f**25f**25#“»85#“«9/ 

data I yoin.l = 1.811/56.3. 56.7. 5S«l.S3.6.$2.ltS0«S|A9.5t86,2t 



126 

127 

126 

129 .. 

13C 

131 

196.8,95*2.93*9,91*1,38*8,36*6,39*5,32*9,31*2.29.1,26*6,29*0,19*5, 
7,3.9.c.Ri.-7.li.-1.7.-?1.5.-2S.2,-2 7*.r’28.7.-3n*5*-32*3*-39*l.-36** 



3-38* ,-90* ,-91*7,-93.5,-45*8 ,-98 *1.-50* 9, -0*0/ 
llATA irrNAtaTl f~l 15I/15*D./ 



data IFCMA31 1 1 , 1 = 1 , 151/15*0*/ 
(ilTk lrC*A3l I 1 . 1 = 1 . 151 /15*u*/ 



132 

133 
13** 
135 .. 
13*> 
137 

OAT* IZ II 1 , 1 = 1 ,91 1 /-35000* , -39oO..“ 3800. ,-3700* ,-3600* ,-3500*0, 
1-3900*0 T-3Tnn.n.-3?OU.O*-3iOD*Q.-3O0Q.Dt-29QQ.Q.-26na*Q_,-27Dfl*Qi. 



■ ~ 2-2600^0i-2500*0, -2900*0, -2300*0, -2200*0, -2100*0, -2000*0, -1900*0, 
3-i6QD.g.-17o{).n.-16nn.n,-lSOO.Q.-1900*0,-13DO*0.-1200»0,-lloo.n. 



9- 1000 *0,-90o7o, - 800*0, -700* 0,-600*0.-500* 0,-900 *0,-300 *0,-200*0 , 
5-|nn. □•□*□/ 



138 

139 
IMF 
IHI 
IN? 
1<*3 

DATA IT inX II 1,1 = 1, 91/0 *0,0*2, 9*2, 9 *9, 5. 0,6*0, 7*0, 8*0, 90*0/ 

HAT* iPFRAl?iTl.Irl.l21/n...7l7*.389.*525*.6a5*l**0**5**10**lS*, 



*25* ,90*/ 

DATA IZCAA2II1. 1 = 1*121/6*0*. *1 19, *19. *135**13. *117, *085/ 



data IZCNA2II1, 1=1, 121/6*0*, 0*. *03, .055, *085, *175, *29/ 

DATA IZCHA2lll.I = 1.121/6*0*, 0 * , - * o35 , - . 065 * - * 1 , - * 2 , - * 28/ 



!•«% 

IMS 

1*16 

1M7 

DATA IZCAA2lIl,I=13,291/*212,*225,.225,.223,*21 ,*16,*289,*3,*325, 

*.325. *3*. 2/ 



DATA |ZCNA2III,I = 13,29 1/0., .o3 . * 07 1 . * 1 1 5 , * 2 , * 2 7 . * 0 . *02 8 , *069 , * 1 1 , 
**185**27/ 



1*IP 

_199 

13F 

151 

152 

153 

DATA l2CMA2in,I = 13.28>/0..-.04,-.08,-.13.-.23,.,3i,0.*-.033f-.075 



DATA (2CAA2flltl = 2S«36)/.3;'7,.3SS,.8l5».8lS».39v,32».55t.SSt«S5t*5 



DATA <2CNA2( I ) , 1 = 25. 36 1 /O. • .02 , .052 , .095 , • 1 75 « • 275 . ,0 . • 01 1 • 032 t *09 



1&9 

|bS 

DATA l2CHA2lil,l:25.361/0*,-*02,-*06,-*115.-*2.-* 295,0*, -.013, -.03 
*7. -.091. -.175. -.28/ 



156 

157 

^^Vo9 l^TTir 

j=i 

H0000990 

H000U950 


158 

159 

FCIl,ji=0*0 
b09 CohTTNUF 

H00009e0 


UD 
lfcl_ . 
167 
163 

NM:0 

AA0=1*3302117l-0Z 



AAi=-8*8502069E-05 

AA2=-9*2193056E-l9 



16*1 

IbS 

AA3:5.9517557E-13 

AA9=-3*9799789E-17 



lb6 

167 

AA5=7.8771273E-22 

TpR=O*0 



168 

169 

Read 3,in,ni,n21 

READ 1 , IRAPP.RETR.ZMRPl ,ZMRP2 ,NPARA I 



1 7f. 

PRINT 1 , IRAPP.RETR ,2MPP1 ,2MRP2,NPaRA 1 






flCAO l«IOT,DCL 

.PRINT l,lOI,QELfLI 

READ l,|Xl,Vl,21,UI,Vl,Wlt 

bead 

READ 1« IX3.Y3,Z3,U3,V3.u3> 

U P 1 1 ai .W l.TnTi.PHll.PSin 

read |. IP2.02.R2.THI2,PH12,PSI2) 

-Pl:PnS7.3 

Q1=Q1/S7.3 

Wl=Rl/57.3 

P2=P2/57.3 

7..,5 

R2=R2/S7.3 

JMIl=THll*O.OllAi 

TH12rTH12*0.017««i 

pHll = PH|l*u.0n5S 

PH12=PHI2*0.o17%N 

PSll=PSll*li .017N«i 

P$I2=PSI2*0.017*MI 

read 1. <P3.«.BJj^l«IlxPHI3xP513» 

READ l.(XlXl,Xiri,Xl21,Sl.TL.XH|> . 
head 1 . ( X I x 20 , X ll?QxX lZi_,S12 ^TPJL , X XN 2 l_ 
READ lftXlX3«XlV3.Xl23»S3.TPL3tXXM3l 




READ 1. (XLX2.XLV2.XLZ2*XLBXl,XEBTl.XL62il 
™RE#0 l,lXHZlTj_XHZ2.Liil!ilTj.XHjixXHJ2»XN33!!_ 
READ 1.I6«RALPH| 

PRINI 3,«N|N1,N2! 

PRINT 1.IX1.T1»21.U1,V1,W1) 

PRINT 1.(X2.Y2.?2.U2.V2.M2» 


PRINT 1.IX3.V3.23.U3.V3.M3I 

PRJNT 1. (PliOliEiilHUjPNIlrPSIll 

PRINT 1, |P2,Q2.R2.TMl2.PHi2.PSI2t 
PRINl.l,«P3f(l3^R3,TH13i£HI3jPSI3l_ 

PRINT l.IXIx'l.xIyl.XlZl.Sl.TL.XNil 
PRINT 1.»XIX2D,XIY2D.XI22 .STZ.YPL .XXW2I 


PRINT i.lXlX3.XlT3.XI23.S3.YP|,3.XXM3l 
_ PRINT 1 , IXLXltXlYl .XL ZltX LAXl.XLAYl.XLAZl I 
PRINT 1 ,|XLX2.XLV2.XLZZ|XLBXl,xL8Tl.xLBZir 
_ PRINT 1 , 1XH21T.X M22T.XH 23T.XM31.XH32.XH33I 
print l.lCl 

_3 format t3l2» 

2 FORMAT ISE15.6I 

1. FORMAT. I6E12.6! 

XMAALlli: .A 

XMAAL 12 1 =_.i» 

XMAALI3I = .9 

_D^5Q_J .=_5jl18 

XMaaLII) = AMI-31 

. N*Z«H = 3 _. 

NAZI2I = 15 

„ MSZin = -1 

NSZ«2» = “I 


PoRO=lS.b 

P0R0=16.Z 
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4 ^ 

4 ^ 
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2ke K0V2 : 1. 

2^9 _ - ... B0V3 :_ 1 . 

23C X2FR = *2 

231 12P9_^_Y2 

232 Z2PR r 22 

_ 2J3 ITTi r 0 

23^ 1TT2 = 0 

235 r4«Z2llt : b 

236 N«22«2l : 6 

237 _.NSZ2|1| = -1 

23P NS22I2) = -1 

239 S Sq - . 1 

29C SSU=0.0 

2«*1 ... .SSO=.Oi. . 

2«? SSI = .7 

2%3 SSi=0.75 

2H<i SS1 = 1.U 

2<*5 »l=Qt8 

296 SS2 = 1. 

297 _ _ m_: 9. 

29fi 711=6.0 

299 771 = 12.0 

250 772 = 20. 

251 Sn107 = .55»ST2 

25? 7EMP0 = .65 ♦ PORO * YPl 

253 CPl=D.O 

259 CP2=o.O 

255 CP3=0*O 

256 CP9=0.(i 

257 CR5= 0.0 

256 CP6=0.0 

259 00 799 1 = 1, 8_ 

26f CPlzCPl.OCOSXlil 

261 . _ CP2=CP2«0C0SYJI1 

262 CP3 = CP3«OCOS2l 1 1 

26! CP9 = CP9»ILCY1H»DC0S2« I >-LC2<l l*OCOSY« I) 1 

269 CPS=CP5*|LC2l 1 19DC0SXI I I -LC X 1 1 1 9DC0S2 1 1 1 ) 

265 CP6 = CP6»|LCX7ll*OCOSYIJi_-LCYmiOCOSxmi 

266 799 C0N7INUE 

26 7 Q11=S1N >PH1 1/2. >»C OStPSH/ 2.>»COS7 7Hll/2. l-SlNlPSIl/2. I*S1M 7HH/2M0001 050 

266 •. l♦c6s^PMIl/2. » ' Hbddl060 

2 69 Q 1 2J S.1 NJ PH12/2. 79C0S7PS12/2 . >*C 0SnHI2/ 2.1- SIN I PS 1 2/2. » *S 1N7 7H12/2M00Q 1070 

277 *. I*C0S(PhI2/2. I N0001080 

271 _ 013=SlN7PH13/2. »*COS7 PS1 3/2. )»CQS7 7H13/2. I-S1N7PSI3/2. >*S1N|7H13/2 M000 1 09Q 

. 272 *, »*COS7PNi3/2. » MOOOilOO 

273 _ Q21=SlN7 7HIl/2.t*COS7PSH/2.l*cOS7PHH/2. ^♦SlN^PSH/2. l♦SlN^PHll/2^i000^110 

279 *. 1*0031 7HU/2. 1 H0001120 

_ 2 ^ 22 = S INI 7 HI 2/ 2. l*C0SlPSI2/2. >»C 0SIPH1 2 / 2. > ♦SINI PSI2 / 2 . I»S1NIPH12/2H0D0H30 

~276 *. I*CoSI7HI2/2. I M0001190 

277 _ 023=S1NI7hi3/2. »*C oS IPSI 3/2 . »*C0S I PHI 3/2 . I ♦SIN I PS 1 3/2 . >*S1 N IPHI 3/2M000USQ 

278 *. I*C0SI7HI3/2.) N0U0n6d 

2 79 _Q3I = SlNIPSH/2. l«COSI7HIi/2 . I *cOS I PhI I /2 . » -S IN I 7H 1 1 /2 . » *Sl N I PH 1 1 /2H00 Q 1 1 7 0 

260 *. |*C0S(PSII/2. I MdOolIBb 

281 _ 0 32 = SINIPSI2/2. I *C0S| 7H12/ 2. I*C0SIPHI2/2. I-SINI 7HI 2/2. »*SINI PHI2/2M0O0II90 

2b2 *. l*C0SIPSi2/2. I ' H0001200 

28 3 __ 0 33 = SIN IPSI3/2. l*C0SI7HIJ/2. I*C0S7PHI 3/2. I-SINI 7HI3/2. >*SIN|PHI3/2H00CI 2 Iq 

269 *, )*COSIPSI3/2. > * - - - Hbi,0i220 
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7aS UAl=C0SIPSll/2 . I*CoSITHI 1/2. >*COSII 

28t .*.I•S 1 N 1 PMI 1 Z^.I 

2B7 GA2=COS(PSl2/2. l•C0S(THI^/2.IAc0&<■ 

288 *.1*SIN|PHI2,/2.| 

289 fl«i3=COSIPS13/2. »*COStTH13/2.l*COSII 

.VC •jJL?SIlMPHIiZ2.tJ 

29 1 Nc0=30 

292 .-_K50=0 

?93 1 : 0.0 

299 11=1 

291 TIME : 0. 

.-vt PELl^Ql 

297 DO 78 I:it3 

298 0078 j:1.3 

299 E3(I.J|:0.0 

300 78 CQNjINuE 

301 13(1.11:1.0 

3U2 LliL.2JL=l*iJ 

3Q3 E3(3.3i:l.u 

309 . XLl II .1 »=0.D__ 

30? *lH1.2»:-XL21 

308 _ _ XU(l,3»^*Lri 

307 XL1I2.1 1?XL21 

308 ll4.!2t2»=JMi 

309 XHl2.3» = -*LXl 

310 XLx( 3,1 »:-XLTl 

311 XL1(3.2>=XLX1 

312 _ XLl|3,3l:0.0 

313 XL2li.l|=0.0 

31? XL2 II.2C -X LZ2 

315 Xt2ll,J»:XLt2 

3U XL2C2,l»:Xc2^ 

317 XL2i2.2>:0.0 

318 XL2(2.3I:-XLX2 

319 XL2(3,1>:-XLV2 

32C X L2I3.2(:XLX2 

321 xl2l3.3):o.O 

322 _ XlA1I1,1I:q.O 

323 XlAHI. 21 : -XLAZl 

329 _ **-Al C 1 , 3i:XLAYl 

325 * ' XLAll2.li=XLA?l 

325 XLAl ( 2, 2 1 : 0.0 

'327 XLAr(2,3i:*XLAXl 

328 _XLA1I3,1»:^XLAy1 

329 xLaI I 3,2):XLAxi 

330 XLAl I 3t3l:U.O 

3ii XLBl(i.l):U.O 

332 XLBH 1.21 : -XLP21 

333 XLB1|1,3|:XLBV1 

339 XLbll2.1):XLB21 

335 xL6li2.2i:6.0 

338 iLBH2f 3»^XLBX1 

337 ~ XLBliS.ic-XLBtl 

338_ xLB1I3. 2I:x lBX1 

339 XLBi(3.3i:u.b 

39C __7t-05 CONTINUE 

391 ■■ XGl(l,li:0.o 


l|PSIl/2. 


l«SlN(pS12/2. 


KPS13/2. 


l«S IN (THll/2flO 001230 

H0001290 

19S1NITHI2/2H00012S0 

Honnl26o 

I9S1NITHI3/2N0001270 

W0001280 

HOOq}29o 

eotoiinQ 

(10001320 

M000 133 0 

HbOQl39o 

(10001360 

PflOCLlAlfl 

h0001380 

H0001390 

(10001900 

H0Q01910 

(10001920 

HQflQJJifl 

(10001990 

HQnDi95Q _ 

(10001960 

H00019T0 

H0001980 

H0g01990 

HOOOiSOO 

HOOO lSl Q 

(10001520 

(10001530 

(10001590 

HQOOlSSO 

(10001560 

W0001 570 

(10001580 
M0001590 

Hddbi6uo 

__ H000 1610 

(10001620 

H00D1 63Q 

' M0doi690 

(10001650 

H0001660 

_ M0 U01 670 

N0U01680 

M0001 69Q 

HOOOiTOO 
H0U01710 

H000lt20 

M0001730 

HeoBi79^ “ 

**0ugi750 

(1000176 0 

M0U01770 

N0001780 

_ _ H0001 79Q 

HOOOISOO 


o» 






0*TE 01U379 PIGE i 

3%r 

313 

XGH2, 17=0.0 
xGl 13.1 )=1H19G 


hoooUTo 

tlODClS20 


311 

315 

XG2I1 .1 7=0.0 
X62(2.17=0.n 


M0001830 

HOOOlalO 


34f 

1«7 

XG2(3«17 = XXW2 
XG3I1.1 7=0.0 


Noooieso 

n0U0186n 


31I» 

319 

XG3l2tl 1=0.0 
XG3(3,i7 = XXy3 


H0001870 

N0001880 


35C 

351 

00 93 1=1,3 
PO 93 J=l,3 


MqoO^BIO 

.0001900 


352 

353 

B13(l,Jl=*e3li,J7 
81111. Jl=«LAlil.J7 


h0U01910 

N8M1920 


351 

355 

615(1, Jl=0.0 
M ,Jl=a.O 


80001930 

H0U01910 


35t 

357 

B23II,J7=0.0 

( I , j i:0.0 


N0001950 

H0OD1960 


358 

359 

B25I1,J)=-E3|I,J7 
B26II .J1 = -«LB1(I.J1 


H0U01970 

8000198 


3bC 

361 

BD13II , j7=0.0 
BD19( 1 . J7=0.0 


80001990 

80002000 


367 

363 

BDi5l I, J7=0.0 
6D16II,J7=D.0 


80002010 

8Q002020 


361 

365 

B023I1,J7=0.0 

GD21(I.J7=U.a 


80002030 

80OQ2010 


36f 

367 

B025II,J7=o.O 

B026(l.Jl=n.0 


80002050 

80002060 


368 

36’ 

93 CCmTInUC 
DO 7 1=1,3 


80002070 

8000208a 


37C 

371 

00 7 J=7,9 
B( l.J7=B13( l.J-67 


80002090 

80002100 


37? 

375 

G0il,Ji:b013(l, J-67 
7 CONTINUE 


80002110 

80002120 


37* 

375 

00 8 1=1,3 
Oq 8 J=lo,12 


80002130 

80002110 


376 

377 

0(l,J7=B11(l,J-97 

B0(1.J7=B011I1.J>97 


80002150 

80002160 



378 

379 

380 

381 
362 
38 3 

8 CONTINUE 
..0 9 1 = 1.3 


80002170 

80002180 


00 9 J=13,15 
B(1.J7=B15I l.J-127 


80002190 

80002200 


bOll,J7=BDlSll, J-127 

9 continue 


80002210 

80002220 


381 

1^0 10 1=1,3 


80002230 


381 

DO 10 J=16.1B 


80002210 


38t 

BI1,J7=B16(I,j-157 


80002250 


387 

Ifi CONTINUE 


80002270 


388 

00 13 1=1,6 


80002280 


389 

DO 13 J=7.9 


80002290 


39C 

b(l,J7:B23ll-3,J-67 


80002300 


391 

bOll , J7=B023(1-3,J>67 


80002310 


392 

13 continue 


80002320 


393 

DO 11 1=1,6 


8000233a 


391 

00 11 J=10,12 


80C02310 


395 

blI.J7=B21ll-3,J-9l 


80002350 


39f 

60I1,J1=B021(1>3,J-97 


80002360 


397 

11 CONTINUE 


80UC2370 



39P 00 15 I:**,*, ■ - . . M0002380 




00 i5 

eilf JI=B2Stl-3 .J-12t 

BDIl tJ>=B025II-3,J-12l 

15_C0hT1nUE: 

00 lb 1=«,6 

.iL=l6.*l£ 

BII,JI = -B2(,(I-3iJ>15t 

BDll,Ji:BD26ari4JTliJ__ 

16 CONTINUE 

DO J9 

00 39 J=l,18 

HIlxJJ =J) f 

39 Continue 

fmtii=xMi 

NI2,2I=XN1 

NI3r3»=XHi 

H(i|,M = XIXl 

IU5x54i>ll II 

MI6f6l=X121 

f1ll3>13i:XN31 

Hllb«16l=XH32 

Ha5.15»:XN33 

N{16»16i=XlX3 

KUlsllllXl T 

NllSf 18IIXI23 
UO S3 1 - 1,19 
00 S3 J=lil8 
HI II ,JI=0> 

53 CONTINUE 

Hill, II = 1 ,U/XH? 

HI I2t2i = 1.0/XH1 

HII3,3I=I.0/xH1 

Mill , 61 = 1 . 0 /xixi 

HiiSi5>=i.o/xi;i 

HI I6t6> = l.u/Xizi 

HHI3.13)=I,0/XH31 

Hljl6, 161=1. 0/XH32 

HI 115,151 = 1. 0/XK33 

HII16. 161 = 1. 0/XIX3 

HllIT, 171=1. 0/XIT3 

HI I 18,l8> = l*b/XI23 

0 SI 1= 1.? 

DO 51 J=l,3 

Hlll,JI=g.U 

h2(I ,J>=0>6 

. _h3«l ,JI=0*0 

1 CONTINUE 

H lii , 1 »=XH1 

H1I2,2I=XH1 

HII3,3i:XHl 

H3ii,n=xH3i 

__H3I2,2)=XH32 

h3I3,3I=xH33 

Da^L 1=1,3 

Do 27 J=l,3 

X jlll,JI=0.0 

X 12 1 1 ,J 1 =0*0 



cx> 
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957 

X13I1 »Ji=0.0 
27 CONTINUE 

H6S6J1SQ 

M0003120 . _ 


958 

959 

X11I1,1I=XIX1 
Xli l7.?l=X171 

r.C'^03130 

N0003190 


9bO 

96L 

»Ill3f3) = )(121 

Xl3(l .1 )=X1X3 

M0003150 

M0003190 


96? 

963 

X13|2,2I=X1V3 

Xl3l3.3i=X123 

M0003200 

M0003210 


969 

965 

FGI9, 11=0.0 
FGI5.1l=o*n 

M0003220 

M0003230 


9 61 
9fel 

F G 16 f 1 1 - 0 *0 

FGIlD.l irn.n 

M0003290 

Moeejsso 


968 

969 . 

FGIlltlt=0.0 
FGI12.1 1=0*0 

M000326U 

M0003270 


9 7C 
971 _ 

FGU6.1 1=0.0 
FGI17.1l=Q.n 

M0005280 

M0U03290 _ 


972 

973 

FGI18.1 1=0*0 

IBC IrtO=9 

M0003300 


979 

•»»5 

976 

977 

GO TO 9 
*>06 1ND = 1 

M0003320 

M0003330 


560 CONTINUE 

CALL RUNGEIIND.TlHE.DCLt.Sc.QI 

M0003390 

Ma0Q3350 


97B 
9 79 

9 00 1 1 =- 1 >R 1*02 l«Oi*() 31 -Pi 9091 1*0.5 
7?=-?1*0.3098 

M0003360 


«BC 

RHO=l 1. 16790729*ExPf AA0«2?*( AA 1«Z2*|AA2»??*I AA3*?Z*IAA9«22*IAAS> 1 > 
nm*o. 00199 


962 

98? 

OElDI = TEMPO * ISOPTtSSlI - SQRT(SSO)| * KOVl 
0ELn2 = TfNPD * (S0RTISS21 - S0RTI551I1 * X0V2 



969 

985 

0EL03 = TEMPO * 1 1 - S0RT1SS2II * N0V3 

OELISO : OELOl * OELOl 



966 

967 

OEL2SO = 0EL02 * DEL02 
0EL3S0 = 0EL03 * DCLD3 



966 

989 

IF (TIME.GE.TT1 ) 60 TO 5007 

XYZPRM = (XZ - XZpRJ*^ ♦ (Y2 -j2PR>9*2 ♦ (Z2 - Z2PRI**2 



99P 

991 

PER = SSU ♦ (SSI - SS0»*XYZPRM/0EL1S0 

SDDOT = (2.*U?*SOTOT*(SS1-SSOI*SQRT(XYZPRMI)/OEl1SO 



992 

993 

IF (PER.LE.SSII GO TO 7020 
PER=SS1 



999 

995 

SOOOT = 0. 
GO TO 7020 



996 

9?I 

5J07 IF (SSl.EO.il GO TO 7020 

IF (TlME.6t.TT2l 60 TO 1070 



996 

999 

IF (ITT l.EO.n GO TO bOOT 
X2PR = Xf 



5LC 

561 

Y2PR = Y2 
22PR = 22 



56? 

363 

ITT! = 1 
6607 CONTINUE 



569 

565 

XYZPRM = 1X2 - K2PRI**2 ♦ (Y? - Y?PRI**2 ♦ (22 - 22PRI**2 
PER = SSI ♦ (SS? - SS1I9XY2PRM/0EL2S0 



566 

5U7 

SOOOT = (2.*U2*S0T0T*|SS2-SS1I*S0RT(XV2PRMII/01L2S0 
IF (pER.LC.SS2l GO TO 7020 



568 

569 

PER=SS2 
SOOOT = 0. 



511’ 

511 

GO TO 7020 

lu76 IF (SS2.EQ.il GO TO 7020 



512 

IF (ITT2.E0.il GO TO 1170 






• 9 


f 
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513 

5»* 

X2PR r *2 
Y2PR = V2 



SJ5 

S16 

22PR : 22 
1TI2 : 1 



517 

51B 

1170 CONTINUE 

XT2PRM = 1X2 - X2MI*«9 ♦ IY2 - T2rtl*«2 ♦ 122 - Z2MIM2 



519 

5£0 

PER : SS2 ♦ (1. - 552 l•XY2PRN/0CL350 
5DDOI = «2.*u2*50T0TV(i.-S52iVS0RTlXY2PRHn/DEL350 



521 

522 

IF <PcR.LC«l«Oi Go ^0 lO^O 
PER=1.0 



523 
. 529 

5DD0T = 0. 
7U20 CONTINUE 



525 

52b 

9905 52:PER*5T2 

S0R96 = PER95DT0T 



527 

52P 

RALPHA:RALPH«YPL*( -0.2 2*PER *0.221 
ZHAPP=RH0v50RA6V5QRTISDRA6l*RAPP2SQRTINPaRAl 



52V 
53C 

XHAPP : <-.8*PER * 1.2I*2MAPP 
YHAPP : xHAPP 



531 

532 

h20oT = 1.5*RHO*5DOOT*SORTISOR9G|*KAPP 
XN21=XHAPP*IXXU2/61 



533 

539 

xn22:tmapp*ixxu2/g> 

XH23=2NAPP* IXXW2/6I 



535 

536 

xix2=XIX2D*YMAPP*RAPPARAPP 

XIT2=XIT20«XNAPP*RAPP*RAPP 



537 

538 

M7t7 | = XN21 
HU,ei = xN22 

h00025Bq 

H0002590 


5 3V 
590 

H<9,9UXH23 

HI10.I0I:XIX2 

H0002600 

H0002610 


591 

592 

Hlllflli:xlV2 

HI12,121=X122 

H0002620 

H0002630 


593 

599 

Mil 7,7I=1.0/XN21 
HUB. 81=1. 0/XH22 

n00028oO 

H0002810 

- 

595 

596 

HI i9,vi=1.0/XH23 
HIIIO. 101 = 1. 0/XIX2 

H0002820 

H0002830 


597 

598 

Hlllltlll=i.0/X1Y2 
HH12. 121=1. 0/9122 

H0002890 

H0002850 


59V 

550 

H2(1.11=XH21 
H2I2.21=XH22 

H0003010 

HP0P3O2O 


551 

552 

H2I3.3)=XH23 

X12(1.11=x1x2 

H0003030 

H0003160 


553 

_55? 

X12I2.2 1=X1Y2 
X12I3.31=X122 

H0003170 

H0g03l80 


555 

556 

0012=-«-R2*022*C2*032-P2*0921*O.S 
QD13 = -1.-»»3*OJ3*OJV033-P3909 3190.5 

H0003370 

H0U05380 


557 

558 „ 

0013=0.0 

0021=-|R1*Q11-P 1*0 31-01*091 1*0. 5 

H0003390 


559 

S6D 

OD22=-IR2*012-P2*032-02*0921*0.5 

U02S=-IR3*ul3-P3*033-03*0931*0.5 

HU00390 

H000391 


561 

562 

0023=0.0 

0031 =- 1 -0 1 *0 1 1 *P 1 *021 -R 1*09 11 *0.5 

H0003920 


563 

569 

0032=-l-02*012*P2*6?2-R2*092i*0.5 

QD33=-I-03*013*P3*023-R3*0931*0.5 

H000I93U 

N0Q03990 


565 

566 

0033=0.0 

Q091=-IP1*q11*Q19021*R1«031 1*0.5 

H0003950 


567 

568 

Oo92=-IP2*012«02*022*R2*032 1*0.5 
OP93=-«P3*013*03*023*R3*0331*0.5 

HOOo39bo 

H0003970 


569 

0095=0.0 





oi 

o 




7u01 CONT 

OJU 

Gill 

0 J U 

G1 12 

UJl 

G1I2 

bl«3 

G1I3 
CH3 
7U02 CONT 


588 

G2|3,2)=: 

589 

02(3, 3I=- 

590 

7U03 continue 

S91 

03(1. 1)=( 

592 

G3(1.2I=: 

. 591 


599 

03(2,11=; 

595 

0312. 2I=> 


1 1 =« rQ2i*^L2_r9ll**2 

2) :2.0*(011*021*031«0mi 

3l=2.0<!(Qil^l^«l*Q!Lll 

ll=2.0*l011*021>031*0in 

2I=-Q 1 !♦* 2 *021**2 -031**2 ♦0«il* *2 

3I=2.o*I 021*0 31*Q1 l*Q<i II 
ll=2.b*IQll*J13Llft21*4'U„» _ 

2l=2.U*IQ2U0 31-011*Q<ill 
3i:-01l**2 -fl2l**2_l_Q3l**2 
NUE 

ll_=_S 12**2 - Q22»*2 - 032 *» 2 ♦ 0«i2 <»2 
2l=2.li*IQ12*022«032*QH2l 

3) =2.0*tl312«Q32-fi22*812l_ 

1 I=2.0*l0l2«022>032*0<i2l 

2 1 :-012**2- tOZZ**!- -xjJ2**2_*QH2**2 

3i:2.0*IQ22*032*012*0H2l 

i.U2itn*lQ12*fl32*tl22*fl‘l2l 

2I=2.0*I022*0 32-Oi2*0<»2> 
3l=-012**2.-022**2_«a32**2 ♦812**2 


!>9E G3I2, 3I=2.o*I02 3*8 33*013«813I 

597 G3I 3. 1 1 =2»0*l 013*8 3 3*82 3*013) 

S9« G3I 3 t2) =2. U* 1 023*033-0 13*0131 

599, G } I 3.3)=-013 **2 -023**2 *03 3**2 *013**2 

(•GC 7U01 CONTINUE 

EGl .. __ CALL MTXMPTIG2fGl,6Jijii_r3, 3) 

bu2 call NTXHPVI63.61,G13«3(-3,3) 

603 ..00 91 l = lt3. 

601 00 9n J=l,3 

605 B1H1.J) = G12II. J) 

606 B21II.J|rG13|I,JI 

607 99 CONTINUE __ 

6UP 00 89 1=1,3 

6L9 00.89 j:l,3_ 

61C 01211, JI=-Gi2ll,Jl 

6)1 .6i3iItJ)=T6L31l,4i 

612 89 CONTINUE 

613 . CALL NTXHPYl012.tXLljBL2i3t3|3) 

611 CALL P1TXHPVIG13,XL2,B22t3t3»3i 

615 _ _ 00 95 1 = 1,3 

616 Oo 95 J=l,3 

61 7 feliU iJ I = HB !2il*j£l 

618 01311, ji:-G13II,J) 

619 9b. CONTINUE 

620 XOMi I 1 , ii=u.O 

6£l _ X0Nlli,2l=-Rl 

622 X0Nlll,3l=0i 

623 AQMI2,1) = R1 

621 XONl 12,21=0.0 

625 XOnl I2,3I=-P1 

62f. XONl I 3, 1 )=-01 
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KAIN 


627 
62P 
629 
6 3C 
631 

_ 

*■33 

639 

635 

636 

637 
636_ 
6 39 
69C 

691 

692 

693 
699 
695 
<«9fc 

697 

698 

699 

65D. _ 

651 

652 

653 
6 59 

655 

656 


657 

658 

659 

66r 

661 


663 

669 

665 

666 

667 

668 


669 

67G 

671 

672 

673 
»79 

675 

676 

677 

678 

679 

6«C_, 

661 

bB2.. 

6b? 


X0Hll3.2t=Pl 

X0(Ul3,l|rQ.Q 

X0H21 1, il=U.O 

X0fi2ll,2l:r»2 

X0H2I 1.3>=Q2 

XOH2t2a>=Rl 

XOH2l2t2)=0.0 

-_X0J1^I2.3I = -P^ 

X0H2I 3. 1 ir>Q2 

JOM2l3,21=PJ{ 

X0H2I3,3)=0.0 

XiMjJ 3 ai=0*D 

X0n3l 1.2 I=>r3 

J!9H3I1,3>=03 

X0H3l2«l l=R3 

_X0M312,?|:0.0 

X0H3I2.3|=-P3 

JUIM 3JL3 l1L= -<LS 

x0h 3I 3f2)=p3 

- _*0»i3l3f3i:0.0 

00 88 lri,3 

00,88 J = l,3 

xbn 2 | I t Ji =-xbn 2 ( It J) 

__ *5113 J 1 1 J » =~xO M3 n. j> 

88 CONTINUE 

„C«U HTXHPVIXOH2,Gi2,iAjlti, 
CALL MTXHPT IG12tXOMl,eB.3t3i 

00 81 1=1, 3_ 

00 81 J=l,3 

B011II|JI=AAIl.JI*BB(ItJ> 


81 CONTINUE 

CALL HTXNPT1X0N3jG_13jAAj 3,3, 
CALL NTXHPri6i3tX0NltBBt3t3, 

00 87 1=1,3 ,, 

do 87 J=l,3 

B02HI,J) = AAn.JMBBlI,J| 


87 CONTINUE 

00 82 1=1, 3__ 

do 82 J=l,3 

_B011II,JI^B0lllIiJ» _____ 
bb2iii,j)=-B02iii,jT 
82 CONTINUE 


CALL NTXMPvlB011,XLl,6D12,3, 
__ N T XMPT » 802 1.X L2, 8022 ,3, 

00 80 1=1,3 

_ -P9_8P J = 1 ,3 __ 

x6n 2( i , JI=-X0M2< 1 , J) 

X_0M3l 1 , J l=-X0H3n, J1 

bdll(l,J7=-B011ll,JI 

b021II,JI=-B021lljJI _ 

8L CONTINUE 

PO 5 1=1,3 

00 5 J=l,3 

B J 1, J1 = B 1 1 e I * J I 

BDII ,J|=B011 II , J) 

5 continue 

00 6 i=l,3 
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N0009050 

H000906Q. 

80009070 

H0oa**Q9iL 

H0U09090 
H0009100 
M0009110 
_ «0009i20 
M0009130 
MOoq 9 1 9jJ 
80009150 
80009160 
80009170 
_ 80 009180 
80009190 
8nn092 00 
80009210 

80009220 

80009230 

800092^ 

80009250 
8Q009260 
80009270 
80009280 
80009290 
_ 80009 300 
80009310 
800Q 9320 
80009 330 ’ 
80009390 


80009350 

80009360 

80009370 

8Q00^3Sg 

80009390 

80009900 


80009910 
80009920 
800099 Sb 
80009990 
H0009950 
80009960 





f>89 

00 8 J=9,6 

. 68S„ 

eUtJi = B12(l>J-J) 

888 

BDII iJl=B012II tJ>3) 

bfcl 

6 CONTlNUL 

888 

00 11 1=9,6 

889 

DO 11 J=l,3 

89C 

BII.JI=B21(I*S,JI 

891 . 

fiDlIiJ)=BD2l(l-3.Jl 

892 

11 CONTINUE 

893 . 

J)0 12 1=9.8 

899 

DO 12 J=9,8 

895 _ _ 

Bil,JI=B22ll-3.J-3l 

(98 

BDli.JI=B022(I-3tJ-3l 

897 

12 COmTINUE 

698 

XVlIl, 11=0.0 

699 , . . 

. XVI(1.7I=-H1 

70C 

XV1U.3UV1 

7UI 

xvi i7mi 

7L2 

XV 112 , 21 = 0.0 

7ii3 

X»li2^3i:-Ul 

7o9 

XV1I3.1 ,=-Vl 

7LI 

XVl(3.2l=ul 

7U8 

XV1I3, 31=0.0 


XV2I1. 11=0.0 

7liP 

X*2«l »2*=-N2 

7U9 

_ .XV2ll.3l=V2 

710 

XV2(2tl (=M2 

711 

XV2I2. 21=0.0 

71? 

XV2|2,31=-U2 

713 

X«2l3.i:=-V2 

719 

XV2I3.21=u2 

71* 

XV2I3. 31=0.0 

718 

XV3I1 ,1 1=0.0 

717 

X*H1,21=-N3 

718 

XV3I1,31=V3 

719 

XV3l2,l l=M3 

72C 

XV3I2, 21=0.0 

721 

XV3I2.3I=-U3 

722 

XV3I3,1 1=-V3 

72 3 

XV3(3,21=U3 

729 

XV3I3, 31=0.0 

Ui 

DO B3 1=1.3 

728 

DO 83 J=1,S 

121 

XOHKl.Jl = -XOHlIl.Ji 

728 

83 continue 

72? 

CALL HTXHPTIH1,X0H1,A^3,3,3I 

73C 

LALL HTXHPTIAA,V1H,F1I,3,3,1I 

?n.- — 

CALL HTXHPTIX0Mi,X11,AA,3,3,3I 

737 

CALL HTXMPy |AA,X0Hll,XLlI.3.3, 11 

733 

00 37 1=1,3 

739 

Do 37 J=l,3 

735 

. xoni(i,ji =.r* 0 MKUJi 

738 

X0M2(l,Jl:-xbH2(l,jl 

111 

37 CONTINUE 

738 

CALL MTXMPt(M2,XOM2,AA,3,3,3» 


739 CALL n 

790 ■ ~t*LL HTXHPvU(0M7,)tl?,«A,3,i, 3i 


H0009fc20 

HfiOOlfelfi 

N000II690 

H0009660 

n0009feT0 

nooo96to 

W0QQ 9»9Q 

M0009700 

BPQQ9 71Q 

HOOq972o 

BB8g*T?fl 

H0009790 

B000«7SD._ . 

H000976Q 

- _ »M 0 ? 11 Q 

M0009780 

HOOO«»790 

M0009800 

_ HQ0Q981O 

M0009820 

hoog.9|so 

N0009*90 

H00098S0 

M0009860 

_ H0009870 

MO 669 B 8 O 

H 00Q 9 890 

H0009900 

M00099I0 

H0009920 

MCI? JO 

moo9 99o 
W0009950 
H0009960 

H000997Q 

N0009980 

_ itoogil??-. 

MOODS 000 

M0005010 

H000S020 

MOOOSQSO 

MOOOSObO 

i*ogg^g7o 

MMOS080 


NOoOSlOO 
HOOPS 1 10 
H000Si20 

i’oggsiso 

NOOOS190 


MCU0S160 

H0U0S21(T 
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C*LL HT)mi>Tta»,x01V21tirL?It3t3il| 

f 2 1 1 1 # 1 1= F 2 1 1 1 1 J 0?IO u 1 1 U ♦ P2 *» Ll _ 

f 2H2»1 »=F2H2*n»VNAPP*R*f’P*»-*0M0ll0,l l♦o2•R2 I 
F2I«3tU=F2l!4,j»«_iHaP?*»»P?*«r?^2*P2-Q2»Q2» 

XL*H1,1|=XL2HI,I l♦»MAPP•R*PP♦^-*0H0^8,ll-U2•R2 
_ X L2.U? 1 1 l=jiL2iLiiJ l♦X^UPP»R>PP»l♦XO HD ^T.n♦^^2» Q2 
1)0 28 i:l,S 

b0.26.J = lt3 - - 

X0M2I I , JI=-X0M2I It J> 

AOP 3 ii,j»=-xo«m.ai_ 

COMTINUC 

_tlkk Ml XHPT<W3. X 0H3.Aa.l.l.3» 

CALL HTxnPyl««.VjM,F3I,3,3,ll 


«U2*P2I 


><0U03230 

H0Q0s;«9 

HOObS2SO 

PQUQUiQ— 

M000S260 


NOOOS260 


call hTXHPT 1X0(13. .. .. 


W0Q5130 . . — - 

Call mtxhptiaa.xoh31.xl3i. 3.3.11 

. - _ DO 25. 1 = 1 1 3 . ... 


H0005390 

N0D05350 


DO 25 J=l,3 

JUMLJi I.J I=“H Pr 3 Uij; — 

CONTINUE 

cxtt NTXHPy (61.161 tieUliLtiL . 
CALL MTXnPVIG2tX62.r62t3.3in 
CALL MTXMPT(63tX63.rK3.3.3.n 
CONTINUE 

CALL >BLI2.WD.Z2.Ni.UM02l 

NEXT 2 CAROS ZERO PUJ_JME NlNOS 

UWD1=0.0 

UND2=0«0 

UHDlrO.O 

MW02=0.0 

hMC3=0.0 

vMonuuol 

'vho3=UWD2 

WIN1I1.1I=UMP1 

MINI (2. pzVWOl 

■JiwM 3 i n=wwoi 

CONTINUE 

CALL MTXMPt ( 6ItWINltMlWlt3.3.1> 

uwonwiNiiitii 


M000SS60 

H0PQS3TQ 

H0005380 

M0005390 

M0005«00 

M0005A10 

MOOOSAZb 

M0U05A30 

NOOOSAAO 


NMtM»0 

M0005<60 

N0005970 

WOOOSA80 

NOOOSSIO 
WOO 5 510 
WOO 5 5 30 
W00S59Q 
W005550 
MOOOSSGO^ 
M0005570 


780 

VN01=M1N1I2,1) 




M000558U 

781 

hWDi:wiNl(3.n 




~Mb0d5590 


jtiNmjii=u«fi2_ 

NIN2(2.1I=«W0Z 

NlNZ(3tH:i>H02 

CONTINUE 


H0005600 

H00056IO 

N0g05620 

n6005630 


78F 

CALL MTXHPTl62.NlNl.WlN2f 3.3.1 1 



I1000S690 



7E7 

UUD2=W1N2(1.1( 



M0b05650 




WMD2=NIN2(2.n 

NW02=wiN2l3tn 

UIN}I|,1I=UU03 

NiN3(2ti)=VUD3 

6IN3I3.1 l=NU03_ 

CONTINUE 

CALL MTXMPT(63.WI N3.wIN3.3 . 3.1 1_ 
UNOirulNni .1 1 

VM03:WlN3(2.n „ 

ww03rUlN3|3.I ) 


M0b056fc0 

M0005670 

W005680 

M000569d 

Ngug57oQ 

N00057I0 

M0005720 

MOU0573U 

H0ugS79j] 

M000S750 
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79B 

199 

(•ur 
*01 _ 
RC? 

809 

0»j5 

Ai>(< 

B(i7 

BrB 

eu?_ 

PIO 

Bil 

81 ? 

013 

P19 

815 , 

RIB 

817 

RIB 

B19 

B?0 

02J__ 

bi2 

B23 

B2R 

B25 

R?6 

8?7_ 

PiB 

829 

B5o 

Bit 

b3? 

.B33 

83* 

B35 
8 36 
P37 
P38 
_I3*> 


PRO 

8*41 

PH? 

6h3 

PHH 

-Ph5_ 

ETl 

_F32 
F Y3 
- FZl 
FZ2 

FZ3 

8H6 

TLl 

Rh7 

TL2 

BHB 

TL3 

8H9 

YMl 

8 511 

YM2 

«*5I_ 

YM3 

852 

YNj 

853 

VNi 

85H 

TN3 


7bl2 CONTINUE 

UHl=ul-uu01 

Uh2=U2-UU02-R«LPH«*02 

yM3:U3-U«o3 

VH1=V1>VWD1 

¥M2=»2- VMD7« R 0LPH»»P2 

VH3=V3-VWD3 

NMlrNl-MWOl 

MH2ru2-|^h02 

«h3=«3-MW03 

¥ ll:lUHl«UHl«VHI*VHl«WHl9WHi 19*0.5 


V31=|UH3«UH3*VH3*VH3»NH3*UH3I9*0.S 

¥S0uND=0.0D37*2l*1138.n_. _ 

NACh1:V11/VSoUNO 

_ ALi:*RcOSlMHl/¥ll.l9ilAi_ 

AL 2 =«RC0S(UH2/V21 1*57.3 

AL1=JAC0S tM H^/VTll*S7.3 

7U13 CONTINUE 

CAUL TBLN0»rl.3tZCNA,XtliAUi!Al,2XA,NS?,CNBll 
CALL TBLNOI-1 .3,2CM«,XHAAL,NAZ,2XA ,NS?,CHA1 i 
- .. call TBLN0l-lt3.2CAA_tiPAALiNA2,,2XA^NS2,CAAll 

CNa1=C«A1*CNA1*I2mRP1/TL» 

C-*LL TBLNDl-1.3.2CNA2.PrBAL2.NA22.2XA2.i.S22. 


C-*LL TBLN Dl-1.3.2CNA2.PrBAL2.NA22.2 XA2.i.S22. 

CALL TBLNOI-l .3.2CMA2.PERAL2iNA22f2XA2.NS22i 
CALL TBLN0i-1^3.iCAAij.PtRAL7,NA?2,2XA2,NS22, 
CNA2rCHA2*CNA2*l2MRP2/yPL) 

CALL ]Bl! ALtfCMP3|*LJiJltC_NA3> _ 

CALL TBL i «L .rCA A3. AL3.N.CAA 31 

CALJ^_TBL I AL irCHA3>AL3fN,CMA3l 

7U19 CONTINUE 

UAl =.0.5*RH0*V11«V1L*SJ 

OBARI rOAl/si 

gA2;n.S*RHQ*w21*¥ 21*S T2 

0BAR2=QA2/ST2 

^*3=0. 5 •» HO 9 y 31 9 V 319S 3 

U VC t=IUHl*UHl«V HI* VHl 1*90.5 

UWG2=lUH29UH2*VH29WHi»**0.5 

UVG3= |UH39UH3*VH3*VH3l9*b.5 

r Xl--QA1*CNA1*1UH1 /U¥6il 

r X2:-0A2*CNA2*IUH2 /uvG2> 

LX 1=_-Q A39cNA3* <Uh3/UV631 





PAIN 


7blb CONTINUE 

C*Ll NTXMPTIBtNl.>l A.6,ieil <» 

CALL HTXMPY «A1A,B,K,6,-18,6» 

Do 501 l=li6 

00 501 J=l,6 

NINVl I , W I rO.tO 

501 CONTINUE 

_ HlNV(l.l)=j.O 

K1NV|2, 21=1.0 

MNVUt J)zl*IL 

KINVIA, 41=1.0 

1 = 1.0 

KINV(6,6)=1.0 

C»U LINSYSJN,6.t.0LLi.Nlll¥i^L_ 

00 701 1=1,6 

DO 701..J=1,U 

t)ll,JI = -Bll,J) 

__1D JL_C QILLUUL 

CALL MTXNPY<B,HiNV,BiB,*18,6,6I 

DO 703 1=1,6- - 

00 703 J=l,18 

D<1,J) = -B| l^Ui __ 

7o3 continue 

ML* MP YIBI6.BD.T2.ie.6.18> 

CALL HTXNPY(BIB,AIA,TI, 18,6, 181 
. . call. TBL niMX,JHR,HNEj_N2, THRU 

FCll,n=THRl*CPi*RETR 

FCI2,1|=THR|«CP2*RE1R 

r C 1 3, n=THRl*CP3*RETR 
FC<4,1>=THR1*CP4»RETR 

rciSti >= thri«cp5*retr 

_ LC«6,1 I=THR1ACP^*REIR 

00 502 i=l,lB 

J = l_ 

Fsir,jT=r ni,ji*FAii,j)4rG|i,Ji 

5P2 CONT INUE 

CALL NTXMPYIT1,pS,fP 1, 18, 18,11 
_CALL MTXMPTn2,X0NjFR_2jJ8,l8,lJ 
bo 503 i = i , ie 

_J=1 

FR ci , ji=FRi ii ,ji*rR2(i , j) 

50 3 CON TINUE 

DO 504 1=1,18 

Tsiil , j7i? I ?TrJT«f R n , Jr« F A « I , j 

_ 504 CONTINUE 

CALL MTXMPY (HI ,rsi ,X0M0, 18, 18, i 

; THE N EX T 3 C ARDS 7EP0 OUT XON D 

X0M0I13,1 1:0.0 

_ XOHDI 14 ,1 1=0.0 

XOnOI 15,11=0.0 

Call mixhpyiaia,kinv,bbi,-i8,6, 
CALL NTXNPY(BB1,B,CC1,18,6,18I 

CALL MTxmPY«CC1, X0M.X0NC,1m>1 h , 

00 9019 i=r,18 

J=1 . . „ 

XONli ,J>=XONli,JI-XOHCIl ,J| 
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912 

913 

9019 Continue/ 
T : TIMr 

M00q»83Q 

915 

DO 3001 i=ltl8 
DO 3001 J-1 la 

mssiknsi 

916 

917 

XOMt 1 1 t J)=o«o 
3uQl CONTTmUE 

M0006880 

MQOOhugO 

9lP 

919 

DO 3002 1=1»3 
00 300? J=1 ,3 

M0006900 

flPOet’lQ. . . 

92P 

9il 

027 

V2i 

XOHT I i , JI=Gl 1 J, I 1 
X0MTil«6,j*61 = G74ilf7> 

M0006920 

M0006930 

X0HT|I«12tJ*12)=63IJ,ll 
3002 CONTINUE 

M0006990 

M0006950 

929 

925 

926 

927 

C*LL MT*MPy «X0MT ,XON,XOMF ,18,18,11 
Xln=XOHFtl .11 

H00069b0 

V10 = X0HM2 ,11 
21D=X0HFi] .fl 


928 

9it9 

X20:X0MF|7 ,il 
v2D-X0HF|S 11 


9 30 
9 31 

22D:X0MfJ 9 ’ll 
X30=X0HFI 13.1) 


93? 

933 

Y30=XOMFI 19,11 
23DZX0HF 1 |5.l 1 


939 

935 

1ND:1N0*1 

IF 1 IND-91 500.500.600 

h0006B9o 

M0006850 

936 

937 

bOo CONylNuE 

(1 1 RC = AR COS IGl 13.3) l«57. 3 


938 

939 

SSS 1=11-161(1,3 1*61 11, 3 m 9*0. 5 
SSS2=(l-IG2ll.3l*G2(1.3m**0.5 


990 

991 

$SS3=I1-IG3I1 ,3I*C3< l,3m**0.S 
TH11=XT«N2I-G1I 1,3).SSS1) 


"992 

9H3 

TH12=XT*N2«-G2I l,3>tSSS2) 
THI3=XT*N2(-G3I 1.3I,SSS31 


999 

995 

PH11=XTAN2(G1I2,3| ,G1|3,3II 
PHI2 = XTAn2IG2I2.3) .6213.31 1 


996 

997 

PH13=XTAN2tG3l2,3l,G3l3,3l) 
PSll = XTAN2IGl(l,2)fGl(l,m 


998 

999 

PS12 = XTAN 2 IG 2 ll, 2 t,G 2 ll,m 
PSI3 = XT AN2I63(1,2) .631 l.m 


95P 

951 

IF IT-TPR) 9006, 7006,7006 
7bG6 PRINT 1000,7 

M0007ie0 

952 

953 

tpr=tpr*delp 

luOu FORNAT llHO.9HTIHE.3X.ElS.kt 

M0007190 

959 

955 

PRINT 6666, DIRC 

6666 FORMAT ( IH . 9H0 IRC . 3X . F 1 S . 6 1 


956 

957 

PRINT 989B,(CMA1,CNA1,CAA1) 
PRINT 9899, ICMA2,CNA2,CAA2I 


958 

959 

989 6 FORMAT 1 IH , 9HCM A 1 , 3X, E 1 5 .6 , 2X , 9HCN A 1 , 3X , E 1 5 . 6 , 2X , 9HC A A 1 , 3X , E 1 5 . 6 1 
9899 format I1H.9HCMA2. 3x,E15.6,2X,9HCNA2. 3X.E15.6,2X,9HCAA2,3X,Ei5.6I 


9bf 

9bl 

PRINT 7050 

7u5c FORMAT 1 T a , 2HF A , T25 ,2MF 1 , ThO , 2hFG t » 55 , 2HF R , T 7 q , 9H00M6 , I 85 ,3H0MG, 

M0007200 

M0007210 

962 

963 

IT 1U0,2HFC,//I 

PRINT 7051 . IFAIl,l),FIII,i),F6II,ll,FR(l,ll,X0M0(l,ll,X0M<I,ll , 

M0007230 

969 

965 

lFCII,ll,I=l,iei 
7u51 format I7E15.8I 

M0UQ7250 

986 

967 

PRINT 1001,IXl,Yl,?l,THIl,PHIl,PSli) 
luCl FORMAT tlH , 

M0007270 

9b8 

1 2HX1,3X,E1S.6,2X,2HY1, 3X , E 1 5 . 6 , 2 X , 2H2 1 , 3X , L 1 5 . 6/2 X , 9hT H 1 1 , 

M0007260 


H«1H 


0«Te 01057? 


PAGE 


IS 


13X E15.6,2X,1HPH11 ,M,EiS.6,?X|iHPSII« 3X,E15.6I 

PRINT i002.>Ul.Vl,W l,Pl.Ql. Pn 

1002 FqRhAT (1H , 

1 2t4Ul. 3X.LlS.6.2X.2 HY1.3x.C15»6.2X«2HWl. 

l.ClS«6t2X»2HQl,3X,E15.6,2X,2HRlt3X,C15.6) 

EJLIN I IQP3^tXLU 


H0007290 

Ngog73oo 

NOd073io 

3X^tl4*iZ2*.?MPl ,3XH00oZ32fl 
N0007330 


_HDflD7J?JL 

3XtElS.6/2X,1HTHI2,N0U071U0 
^ h0007ll0 




9»? 


lu03 rORHAT ilH 

1- _ _ 3HAL1,3X»I 

PRINT 6001t <X2tV2f 22|THl2,PHI2yPSI2) 

-6oCLI_£0RMAT .UH 

1 2HX2t3X,ElS.b,2X,2HY2. 3X ,C 1 5. 6 t 2X , 2H;2 , 

13X,E t S.6.?X. iHPHI2t3X,ElS»bt2X.1HP$I2,3 X .E 1S,6) 

PRINT 6002, (U2,V2.U2,P2,02,R2I 

6002. format IlH , 

1 2Hu2,3x tClS.6,2X,2HV2, 3X E’1S.6,2X,2HU2, 

l.ClS.6.2X«2H02i3X.Liix6^,2Xj2HR2L3xi£l$«6) 

PRINT 6003,(AL2) 

JMIB3 rORHJ>J„llH_i 

1 3HAL2,3X,E1S.6I 

PRINT. 1010. IRHOjOBARl^aaiRlJ 

1010 rORHAT I1H,3HRH0,3X,E1S.6,2 x«5HQBAR2,3X,E15.6,2X,5H0BAR1,3X,E1S.6I 

HXB = XIX11P1 

HVB=XIV1101 
JL21=*.Ut*Pl. 


m000735o 


MOqoT120 

M00071 30 

3X,E15.6/2X ,2HP2,3XM0007110 
NM076 50 

H000717D 


HXI:Gi I 1, 1 )iHXB«C1 12, 1 I*HYB*G1 ( 3,1 I*H2B 

Hyi=Glil,2)*Hxp«Gll2»2J*MYB«Gl!3,2l*M2B 

H21=G1( 1,31 wHAB* 61 (2,31 •HYB«G1( 3,31 *H2B 

I*H=5QRTIHXIiHXI*nyi*MiJiM2JiH21| _ 

PRINT 9301 , ITAM ,HXI ,HY1 ,HZ1 I 

9 301 FORMAT ( IH , 3HT AM , 3X ,E 1 5 ,6 , 2X , 3hHX 1 , 3X , E 1 5 . 6 , 2X , 3 HHY I «3X«CIS»6 , 2X , 


13HH2I ,3X,E15.6I 
6169 FORMAT I6E12.61 
illi FORMAT I1H,2H$2, 

PRINI..llll,S2 

Nm:nm«1 

PTIMECNMI=TIME 


3X,E15.6> 


PZUNMirZi 

^R7(NM|=FR(7,P 

PFReiNMI:FR(8,li 

PFR9INMI=rRI9,ll 

PSZINMCSZ 

PQBAB1(NMI=QBAR1 


P0eAR2INM)=0BAR2 
_ePl|NM)rPl*51,J^__ 
POl INHi:01i57.3 

PRl (NM)=R1*S7^3 

PZ2 INMirZZ 

PDXl (N M » = X0MF»1,1I 


P0V1(NH)=X0MF(2,1) 

r021«NM|:X0MF(3,1L 

PAL1(NM):a(.1 

PXL2INM I = 

lMTiMe-2O.6>i0O6t 9001 ,9001 

_9oO_6_J1: 1 1 ♦! 

GO TO ZOO 

9u0l CONTINUE 

CALL lOENT (979MCUIPAPEB) 


M0U07510 


h0U07S2u 



C7I 

OD 


P«IN 


DATE 04q379 PACE 19 

1G26 
1027 . 

CALL OUI 83 V l-l,35tTlMELtZlLt-NH,PTiME|P71l 
CALL 0U1K39 1 -1 . 35 . T IMEL . PF RL 7 . -NM . PT IHE . PER 7 > 


ipie 

1C29 

CALL 0mn3W 1-1 ,35,TlHEL,PFRL8,-NP,PllME,PFRj» 
CALL QU1H3V « -1 . 35 , T IHEL .PF RL9 . -NM . PT I ME ,PFR9 1 


1030 

1031 

call 0U183V l-l,35^TlHELtPS?L ,-NM ,PTIME ,PS2I 
CiLL 0U1K3V l-l.iS.TlHELiPDBAlL.-NM.PTIHE.PQBARU 


10 37 
1033 

CALL 0UIK3V t-1.35,TlREL.P0BA2L,-NH,PTlHCfP0B*R2> 
CALL 0UU3W 1-1 .35.TIMEL.P22Lt-NHtPTlME.P72» 


1C3« 

1035 

CALL 0U183V 1-1 ,35tTlMELtPPlLt-NM,PTIMEtPPH 
. . CALL 0UIA3V J-l,35*llMtW^PQlL,-MNtPTl«E,P0l|_ _ 


1031 

1017 

call 0UIK3V t-l,35,TlPCL,PRlLt-NH,PTlMCyPRl) 
Call QU183V 1 -1 . 35 . TImEL .PDX IL • 'NM .PTIHE • PD l 1 1 


1038 

1039 

call 0UIK3V l-l,35tTlPELfPD71L,-NM,PTlME.pDYl» 
CALL 0UIK3V 1 -1 . 35 . T IHEL .PD2 IL • -NH tPT IHE t PnZ 1 I 


io**o 

ICJHl 

CALL QU1K3V |-l,35,TIHCL,PALlL.-NHtPTIHCtPALH 
. . CALL. 0ulR3« l-1.3itLUlEL»PAL2L»rNH*PTlHE,PAL21- 


104? 

1043 

CALL ENDjOB 
CALL EXIT 


1044 

1045 

9101 CONTINUE 

Mnno754o 




dHUC.P RUhCe 


dPRT.S kUNGt 

ruRPUR 27R3DL33 SL73R1 0*l/03/79 18t?9t32 



RUNGE 


31 

3? 

33 

5 * 

35 

Sh 

37 

311 

3» 

<•0 

11 „ 

1 ? 

<3 . 

_H5 

in 

19 

ir- 

51 

57 

53 


7K7lj6*TPr « lU I .RUNGE. 

- ^UBROullNL RUNG E IKUTTA.TlHE, DT »BWtliB>»l 

C RUTT» iS A COnTRqL INTEGER 

C JIML IS Tint or INTtcRAliON 

C DT IS THE TINE INCREHENT 

C — ttJ> AiL_= NU MB ER OF V AR IABLES TO BE INTEGRiTtp ONCE 

C NUVAR: NUMBER OF VARIABLES TO BE INTEGRATED TMICE 

DIMENSION Cll GDI, C^l GDl.Cnil ^ 

1 CD3( 60|,C01( 601, SX( 60I,SXDI 6ol 

iOMMON/FuNI/ xl GDI 

COMHON/FIRT/ XOI GqI 

COMMON/ SE cq/XodUOI 

GO TO (lt7t3t1l*KUTTA 

100-10 l=l,NVAB 

SX II I r X (1 I 

cmi : *01H?JJT 

10 xij) = sx(ii*o.s«ciin 

IFtNQVAR.r O.OlGO TO SO 

DO 100 I=1,N0VAR 

SXOm.:: X0111___ 

conn = XODimOT 

100 xDin = sxoin*o.s*CD iiii 

5u CONTINUE 

UML_^ J IML ♦ D .^>DT 

RETURN 

i 00 20 I=liNVAR 

C2iii = xoin*oT 

20 X|ll_= SXim0.5»C 2III 

IFINOVAR.EO.OlGO TO SI 

DO 200 UltNDVAR 

CD2in = XDOin*OT 

__tOO XDI1»_= SX0lin0.5»C D2in 

Si CONTINUE 

_ return 

3 00 3o i=l,NVAR 

C3tn = X om*DT 

30 xin r sxin*C3in 

IF (NOVAR.EO.OIGO TO 52 

bo 300 I=1,N0WAR 

co3in = xpoii noT 

30L xoin = sxDin*co3in 

5^ CONT I N U E 

TIME = T1ME*0.5*DT 

RETURN 

~1 bo 10 I=1,NVAR 

CIIII =_xpiiIiDT__ 

10 xii) : sxin4iciin*ciin*2.*ic2iiT*c3i?ni/6.~o 

mND VAR,E0,0iG0 TO 53 

00 ■•00 niiNOVAR 

cpiji L_=_*Pbn •♦oi 

100 xDin = sxDin«icDiin*CDiin«2.o*ic52ir)4C03ii »i 

CONTINUE. 

RETURN 

END 


CJ1 

CD 


.<H06,P 





MURPHYBlN?Lf>*1Prs(UI.TBL 


SUBROUIlNt. TbU X.T.XI.N.YI I 

THE DATA MUST 8E IN ASCENDING ORocP 

OltlENSlON Xtll.TIl l 

DO 93 1=1. N 

lflXiJlrXlJ93tl5^9j 

91 0X=Xl-Xll-ll 

D*l=x m-XII-lL 

l,V=TliJ-Y«i-ll 

W=TilTlM.DY*lQlZDAU - 

CO TO 9A 


GO TO 99 

?3 continue .. 

99 CONTINUE 

RETURN 

END 
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•10L»iD 


DATE UA0S79 


r«GE 


NURrHTBIN?l'fi*TPF 

1 

7 C 

3 

A C 

-- _5l JL 


tIUl.TBLNO 

N-J)lMN!iIO.NAL TAfiJLE LO OH-UP R OUTINE 


■j 3/gQ/T3 


IBLNOOlO 

TBLN0020 


IBlnOqSo 

TBLN0060 


7 

« 

9 
1C 
) ’ 

12 
17 
lA 
15 
It, 
17_ 
IB 
19 
2t 
21 
22 
23 . 
29 
^5 
2t 
<7 
2B 

29. 

3C 

31 

32 

33 
39 
35. 

36 

37 
3B 
39 
9C 

J»l_ 


C 

c - 
c 

c 

c 

-L-_ 


NEX7R 


extrapolation control CINiEGErI 

IE_JER(Lt.-MQ_LllBJlPjOLAll5rLJRU« TERkImAT EDJ 
IP POSITIVE. EXTRAPOLATION OCCURS 
I F KEGA T iyL. NO EHIRAPPLAUQN mUH. COllimULPx 


TBLN0070 

TBLN0090 

JLBUNguaiL 


TBLNQlla 

IBI N UIZO. 


c 

c 

c 

G 

C 

x_ 


nD 


table LOOK-UP OINENSION, LIMITED TO 6 CINTEGER] TBLN0130 

EoNE-G8£iilfi_lHiLB KUftBER.Or INDEPENDENT. PARll]E.TLeU_- IBLNOIbCL 

IP positive, independent VECTOR INDICES ARC FOUND TBLNqIsO 

lf-.»l£ilA4HlE^-SAvEQ- VALUES OF INOEEEililENT ItECIllH IBLNOIWL. 

INDICES ARE USED TBLN0170 

..■IAaLE-flf.. DE£LMaLNI PiRaHeUR S CRUl AHBA.TJ TBLND18Q 


TABLE OF INDEPENDENT VECTORS tRC«L aPRAV] 

EEAC>1JEECXQR_HUS1_M-1N .ASCENDING OROER3 

Example of x and ? tablcsc 

- X.-1JB.LE _Z-I#BLE 

A1 CI....C9 FOR B1 FOR A1 

Li 1 » B2 a a 


TBLN0190 

_lfltNO|00_ 

T8LN02I0 

TBLND220 


TBLN0230 

.TBlNg29g 


Bl 

.92-. 

P3 

_C1_ 


a b 3 a a 
a _pi a _a 2L. 
a B2 a a 
a B3 .a..i_ 


C2 


TBLN0250 

TBLNOZfcO 

TBLN0270 

TBLNQ2B0 

TBLN0290 

TBLN0300 


NA 

C9 

LENGTH pF LA&tL INDEPENDENT VECTOR EINTCGEN ARRAY! 

TBLN03I0 

IBLN0320 




X A 

EXAMPLEC NAL3! > 2. 3. 9 
INDEPENDENT PARAMETERS CRCAL ARRAY! 

TBLN0330 

TBLN0390 




NS 

SAVED VALUES OF INDCPCNOCNT VECTOR INDICES EINTEGERS! 
C SHOULD BE INITIALLY SET TO -1 SO THAT INDEPENDENT 

TBLN0350 

TBLN0360 





92 

93 

99 

95, 

9b 

•»7. 

98 

99 
sc 

51 

52 

53 

59 

55 

5t 


VECTORS are tested FOR ASCENDING ORDER! 


ERROR CONOlTlONSt EIF AN ERROR OCCURSf A MESSAGE. IS PRINTED 

AND THE RUN IS TERMINATED!. 

1 . _._A T _t EAST ONE OF THE INDEP ENDENT VECTORS IS NOT IN 

ASCENDING ORDER. 

2. WITH nEXTR set TO ZERO. AT LEAST ONE.OF. IHE_INOEPENi»LNL 


TBLN0370 

TBLN0380 

TBUn0390 

IBLN0 900 

TBLN0910 

TBLNU920 


parameters is out of range of the INOCPENOCNT VECTORS. 
3, TABLE LOOK-UP DIME NSION. Nq , IS GREATER TH AN 6. 


T6LN0930 

TBLN Q990 

TBLN09S0 

IBLNOJW. 

TBLN097Q 

TBLN0960 


DIMENSION 


2(11. x(l). NAIl). XAlII, NS(l). 
R»T1CT5 L,_NGr0UP(.51 . .ITOTISI, WJI321. 


0990 

JileO. 


TBLNOSIO 

TBLN0560 

TBLNObOO' 

TBLN0610 


1 IF (NO .LT. 01 
IF (ND .LT. 71 

2 ttRITC (6.9U2I 


GO TO 60 
CO TO. 10. 


0620 

U690 


1BLND 


1 MRITE l6t9Um 

ilOP - 

C 

10_L1. = .J . ■■ 

tr = NO - 1 

Bn. J5 _J =1 iLt 

L2 = LI « N«(I> _ 1 

If-INSm _.GT. .-Ll_6JB_15_2( 

nSIII = 0 

Jp =_L2_^ 1 

DO 12 J=LliJF 

11 WRITE (6«902l 

MRITt |fc,911|_J^_J 

GO TO 90 

12 CONTINUE 


IXAill .GT. XILl) .«N0. XAIll .LT. XIL2II 

IX«m .GEt laZJLL NSlIJ-=-t2 

IXAdl .LE. XiLlM NSm = LI . i 


1£ JNEXTr) ._33j31LtJ2 
IF INSIII .EO. l2> 

L.= NLLU 

LL = 1 

ira.LT. Lll _GQ_I 
IF ix«m .LE. XILII 
JF ixxm. ,GT. X tH 
L : L ♦ lI 

N Sm = L 

IF IL .GE. LI .«NO. 

60_1Q 21 

WRITE I6t902t 

WRITE I6t930l I 

GO TO 90 

IF I NSm . GE. L?» 

IF INSli) .LT. Ill 


NSin : L2 - 


Ge_iJiL_2:* „ _ 

ILII LL = -1 

in .»W0 . *J«_H_.Lt. XIL^lll 


GO TO 2W 


NSm = L2 - 
NSIll : LI 


H« = NSdl 

KB = NSd > ♦ 1 

IF (NSd I .LT. Lll K« : LI 

IF (NSm .EO. L2» KB = L2 

IF (XIKAI .EQ. XIKBI) GO TO 35 

RATlOin =_iXA|l» - XIK AH / (XIKBI - XIKAII 

GO TO 37 

q.g 

NSd I = K« 

-Ll^E2_±. 1 

CONTINUE 

NGROUPd » ='Nrd» 

NSUM = NAIl) _ _ 

00 WW i:?,LF ■ 

NGR OUPI 1 I = NSm - NSUN 

NSUH : NSUN * NAin 

CONTINUE _ 

ifoTiLFI : i 


OATE 0W0379 


page 


2 




05 


ibiND 


0«TC 040ST9 


P«6C 


11« 

US 

C 

_ 50 00 52- IrZ.lF _ 


TiLNUio ' 

1230 

lib 

111 

J=LF-141 

- llfliJJLr lIQUJtU 


1250 

lls 
liv. 

52 CONTINUE 
C 


1260 

1BLN1270 

uo 

O 

50 

2 0-’ 

6U NF : ?4*LF 
UN : -2 


1250 

1290 

bl 00 59 1 = 1, KF ,2 

IFIRST z 1 


1300 

1310 

HM = NU 4 2 

52 DO 57 J=1,LF 


1320 

fi c ^ 

NH z ?**|j-l) 
MMM ; HM 


1340 

ss* 
^ 1=} ' 
O Ell r 
R 3i 

MWM : HW 

CALL BlNl«HHiD*35»ttjllWMtU - . 



1F«MMN.EQ.0» GO TO 55 
. H6R0UPU1 *1 


1350 

« r 

GO TO bb 

55 i(10N,= NGROUPlJT. 


1370 

13B0 

5 * 

55 IFIRST : IFIRST ♦ « TMON - 1 1 * I TO T I J » 
57 CONTINUE 


1390 

1400 


ISEC : IFIRST ♦ ITOTtll 
NJJJJ : 2lIFinST) 


1410 

1420 

p 

139 

hJ(I«l 1 = 2(ISCCI 
59 CONTINUE _ 


1430 

144fi 

l*iC 

IM 

C 

7C DO 79 1 = 1. LF 


TBLN14S0 

1450 

m? 

143 

KF : 8F/2 
7L-.liG_17_ JzUKF 


1470 

14S0 

144 

145 

NJIJI = mJ( 2*J-1I 
IF «RATI0«ll_,EQ4_0.Qh 0.0. T0J7 


149D 

1500 

14t 

147 

bJIJI : WJ(JI * IWJI24JT - WJI24J-in • RATIOIII 
77 CONTINUE 


1510 

1520 

148 

J49 

79 CONTINUE 

. . 5_ 


1530 

TBLNI540 

Ibt 

lU 

8C ANS = UJIl 1 

RE TURN 


1560 

lb? 

153 

r 

90 CONTINUE 


TBLN1570 

154 

- 11-5. 

RETURN 

_C 


1500 

TBLNlblO 

ISb 

157 

902 format I iHi/iox , 'ERROR OCCURRED IN TBINO ROUTIhE'/lH 
904 .EORmAT UOk.'PIPENSJON QF TABLE EXCEEDS b IM = ',12 

1 

1520 

1630 

158 

154 

911 FORMAT 110* , 'INDEPENDENT VECTOR NO. ‘.I?.' IS NOT IN 
*ER AT. INDEX NO. .',131 

ASCENDING 

ORD 1640 

1650 

IbP 

161. 

930 FORMAT • 1 OX , ' 1 nOE pENDEN T PARAMETER NO. '.f?,' IS OUT 
•HE CORRESPONDING VErTOR IK = Ol’l 

OF RANGE 

Or T 1550 

1670 

lb? 

Ib3 

lOuOO CONTINUE 

END _ . 


16S0 





bHUG,r 

iTan? 



bPRI,5 

XTAN2 



F UMPUN 

?7F30i33 SL73R1 04/03/79 18:29:35 


■ 






1 1 AN?' 


HUNPHTbIM?Cf*TPrsiL) *XT«N2 


.. . 1 . FUNCTION XTAN IU.B t 

2 IF(B .iO. 0,o> CO TO 100 

3 - - ATAN2=ATAN2|A,BI«SItJ 

A return 

i 

i 10 XTAN2r>V0.0 

. 1 —BtIURN 

P ]6 XTAN2=0.0 

9... _ RETURN 

10 IS XTAN2=90.0 

» I auiiRJi 

1? END 


AFIN 



DATt 0A03I9 


PACE 


1 
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